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Abstract 

The Clearwater Formation in the Cold Lake oil sands area of east-central Alberta 
contains the most continuous reserves with the highest hydrocarbon saturation in the Cold 
Lake oil sands area of east-central Alberta. The Clearwater Formation in the detailed study 
area, located at the site of the British Petroleum Pilot Plant (Township 66, Range 5W4), has 
an average thickness of 40 m (120 ft) and is overlain by approximately 430 m (1400 ft) of 
overburden. 

The Clearwater Formation can be informally divided into ten stratigraphic units. Six 
facies associations have been recognized in the interbedded stratigraphic units: heavily 
bioturbated silty shale; heavily bioturbated silty shale to silt and sand; moderately to heavily 
bioturbated sand with silt and silty shale; sand, silt and silty shale; graded sand, silt and 
shale with Bouma sequences; and sand and carbonaceous sand. Two coarsening-upward 
sequences are present in the Clearwater Formation. The facies associations represent 
offshore and nearshore sedimentation in a shallow epeiric sea which received detritus 
from a source situated to the southeast of the study area. Poorly developed very 
fine-grained sand bars and ridges were constructed on top of a muddy shelf and are the 
best potential reservoirs in this area. 

The very fine-grained sands are feldspathic litharenites to litharenites. The 
sediments consist of a complex suite of minerals, including quartz, chert, feldspar, 
igneous rock fragments, metamorphic rock fragments, sedimentary rock fragments, 
volcanic rock fragments, clastic carbonates and glauconite. A dual provenance from the 
Cordillera to the west and the Shield area to the northeast is indicated. The clay mineralogy 
is also complex and consists of varying amounts of kaolinite, illite and smectite with minor 
chlorite. 

The study of the relationship between geophysical well logs and core data 
indicates that correlation of these parameters is difficult. The corrected sonic porosity is 
the best indicator of total porosity and correlates best with the core porosity which also 
detects total porosity. The core and sonic porosity values are overly-optimistic because 
clay-bound water is detected as porosity. 

A shaly sand cross-plot of corrected bulk density versus sandstone neutron 


porosity gives the most realistic values for volume of shale and effective porosity. In the 
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least shaly sands, effective porosity is approximately 12 porosity units lower than the 
porosity indicated by core analysis. The cleanest sands contain 25 percent shale in 
dispersed, structural and laminated forms. Effective porosity is as high as 31 percent. 


Only a small portion of stratigraphic unit F is producibie in well 10-17-66-5W4. 
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I. Introduction 

Due to their economic potential, the oil sands and heavy oil deposits of Alberta and 
western Saskatchewan have received increased attention over the past decade. 

Until recently, most studies have concentrated on the Athabasca oil sands deposit 
because a small part of it can be mined using surface mining methods and the oil can be 
extracted using proven techniques. However, the oil in the vast majority of the Athabasca 
deposit and in all other oil sands deposits is found under hundreds of metres of 
overburden and must be extracted using in situ methods. 

In the case of the present study of the Clearwater Formation in the Cold Lake oil 
sands deposit, in situ extraction requires thie use of subsurface techniques based on the 
study of core and the interpretation of geophysical well logs. With the exception of a 
mineralogic study by Putnam and Pedskalny (1983), no detailed mineralogic study has been 
carried out on the Cold Lake Clearwater Formation deposit. To date, there has been no 
detailed work done on the relationship of lithology and mineralogy to geophysical well log 
response. 

An insight into Clearwater Formation sedimentology and petrology is needed in 
order to satisfactorily deal with a one-dimensional subsurface view of the deposit. The 
vertically, and to a lesser extent laterally, discontinuous reservoir sands, as well as the 
complex mineralogy, reflect intricate facies associations which must be understood so 


that oil extraction techniques can be optimally designed and applied. 


A. Purposes of Study 

One purpose of this study is to define the facies associations in the Clearwater 
Formation by conducting a site-specific analysis of the deposit. From the lateral and 
vertical relationships of facies associations, an environmental reconstruction can be made. 
The patterns discerned from the examination of lithologies and the geoptiysical well log 
response will then be applied on a regional scale. Comparisons with sedimentary models 
can be applied to the deposit. 

The second purpose of this study is to define the relationship between geophysical 
log response and the lithology of the shaly sand of the Clearwater Formation. In cases 


where the Clearwater Formation has been cored, lithologic analysis is relatively simple. 
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However, when wells are not cored or core is missing, one must rely only on geophysical 
logs; and, the analysis is much more complex. 

The geophysical well logging tools respond to the bound water in hydrous 
phyllosilicates, such as clays, mica and glauconite, as effective porosity. This results in 
anomolously high porosity readings. In order to minimize the need for logging core, 


predictive modelling for the Clearwater Formation in this area is imperative. 


B. Location and Geological Setting of Study Area 

The Cold Lake oil sands area in east-central Alberta is located from Townships 52 

to 70, Ranges 1 to 10 west of the 4th meridian (Figure 1). It contains heavy oil reservoirs 

| in all formations of the Lower Cretaceous Mannville Group. In the northern part of the 
deposit, the best, most laterally continuous reservoir is found in the Clearwater Formation. 

The site of the present study is the British Petroleum Marguerite Lake Pilot Plant 
area and vicinity, located in Township 66, Range 5W4M (Figure 2). The Clearwater 
Formation here has an average thickness of 40 m (120 ft) and is overlain by approximately 
430 m (1400 ft) of overburden. Five non-confidential weils with useable core are located 
in this area: 6-3-66-5W4, 6-5-66-5W4, 10-7-66-5W4, 10-17-66-5W4, and 
2-19-66-5W4. The core from these wells was examined, sampled and photographed at 


the Energy Resources Conservation Board Core Laboratory located in Calgary, Alberta. 


C. Methods of Investigation 

Each core was examined bed-by-bed in conjunction with an analysis of the 
geophysical logs. The parameters described were lithology, texture, sedimentary 
structures and macrofossil content. A photographic record was made of each core. 

This formed the sedimentologic basis for stratigraphic correlations and facies 
association analysis. Stratigraphic units and facies associations were defined and 
correlated between these wells. Then a regional study of these facies associations in 54 
wells was undertaken on the basis of we!! log evaluation. Ilsopach and structural maps 
were generated for the various stratigraphic units and the underiying pre-Cretaceous 


unconformity to discern regional trends in a ten township area. 
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FIGURE 1. Location of study area. 
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FIGURE 2. Location of detailed study: structural 
cross-sections indicated by lines A-B 
and A-C. 
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After the initial examination, one well (10-—17-66-5W4) was chosen for detailed 
study. Factors such as completeness of core, presence of representative lithologies, 
good geophysical log control and location near the centre of the B.P. pilot plant area were 
used in choosing this core. 

A total of 43 lithologic samples of 10-17-66-5W4 core were taken. The 
location of samples is given in Figure 3. At least one lithologic sample was taken at the 
top, bottom, and within each major lithology. Maximum spacing of samples was 2.5 m (8 
ft), with an average spacing over the entire length of the core of 0.85 m (2.8 ft). Where 
sand beds were separated by thin silt and shale partings, separate samples were taken for 


oil sands and micropaleontologic analysis. 


Geophysical Well Log Analysis 
All wells were studied in conjunction with an examination of geophysical well logs. 
Logs available for each well are as follows: 


6-3-66-5W4: Dual Induction Laterolog 
Compensated Neutron Formation Density 


6-5-66-5W4: Induction Electrical Log 
Compensated Neutron Formation Density 


10-7—66-5W4: Dual Induction Laterolog 
Compensated Neutron Formation Density 
Borehole Compensated Sonic Log 
10-—17-66-5W4: Dual Induction Laterolog 
Compensated Neutron Formation Density 
Borehole Compensated Sonic Log 


2-—19-66-5W4: Induction Electrical Log 
Compensated Neutron Formation Density 


The geophysical logs were correlated with the core and were used to detect core that had 
expanded after the release of overburden pressure. In the detailed study of well 
10-17-66-5W4, all logs were digitized using a Tektronix graphic terminal. However, the 
digitizing procedure gave readings that were up to 0.5 m (1.5 ft) away from the actual 
depth indicated by the geophysical well log. These inaccuracies were a result of limitations 
of scale on the Tektronix graphic terminal and variations in the geophysical well logs. 
Readings were made every six inches and at every sample point manually. The data were 
analyzed using shaly sand cross-plotting techniques with the scattergram computer 


program developed by Nie et al (1975). The geophysical well log parameters were then 
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FIGURE.3...Location.of. oil] sands..shale (S)...trace.fossil..(7) 
and macropaleontologic (F) samples. 


related to core analysis parameters. 


Laboratory Analyses 

Thirty-four of the 43 samples taken from well 10-17-66-5W4 core were 
submitted for oil sands analysis. The analysis consisted of a calculation of weight 
percentages of oil and water saturation by the use of the Dean-Stark separation technique 
(Kidston, personal communication). The method involves the extraction of the bitumen by 
refluxing hot toluene, measuring the amount of water extracted through condensation and 
obtaining the weight percent bitumen by calculating the difference between the weight of 
the original sample and the dry sample plus water (Eade, 1975). The grain size distribution 
was Studied through both dry and wet sieving. Dry sieving was performed using a series 
of sieves ranging from -4.0 phi to +11.0 phi, graduated in 0.5 phi intervals. The mean, 
median, standard deviation (sorting), skewness and kurtosis were determined graphically 
(Folk, 1968). Histograms and cumulative curves were generated for each sample. 

Thirty-four samples were thin-sectioned, impregnated with blue epoxy and stained 
with sodium cobaltinitrite for the identification of feldspars. Along with seven thin 
sections on loan from British Petroleum, they formed the basis of a petrographic study of 
the Clearwater Formation, in which 300 point counts were made of each slide. 

X-ray diffraction analyses were conducted on the less than 2 micron smear slides 
of 36 samples to study clay mineralogy and to obtain relative percentages of clay types. 

Fourteen samples, 3 from 6—3-66-5W4, 6 from 10-17-66-5W4 and 5 from 
2—19-66-5W4 were submitted for micropaleontologic study by Dr. C. Singh. 
Dinoflaggelate assemblages were studied to discern the geologic age of desposition and 
to aid in the analysis of paleoenvironmental conditions at the time of Clearwater 
deposition. 

Macropaleontologic identification was made of the one complete shell fragment 
found and one sample was taken for trace fossil analysis. Other trace fossils were 
identified through core analysis. In the heavily bioturbated deposits, trace fossil 
identification was impossible. 

A series of parameters were obtained using the equations listed in Appendix lil. 


These parameters were then uséd in the shaly sand analysis. Cross-plots were used to 
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examine and demonstrate the complex relationship between petrography, core analysis 


and geophysical well log response. 
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ll. Stratigraphy 


A. Previous Work 

The "Clearwater shale” was first recognized in the Athabasca area by McConnell 
(1891) and was later elevated to formational status (McLearn, 1917). The Clearwater 
Formation was placed in the Lower Cretaceous Group (McLearn, 1918) along with the 
McMurray and Grand Rapids Formations. Subsequently, the Lower Cretaceous Group was 
correlated throughout Alberta (McLearn, 1931 and 1932). 

The Lower Cretaceous Group in the Vermilion—Wainwright area was divided into 
six members on the basis of lithology and named the Mannville Formation (Nauss, 1945). In 
the Lloydminster area, which is in essentially the same geographic location as 
Vermilion—Wainwright, the Mannville Formation was divided into three units (Wickenden, 
1948). The middle unit is roughly correlative with the Clearwater Formation in the 
Athabasca area and with Nauss’s (1945) Cummings, Islay, Tovell and part of the Borradaile 
Members. 

The Mannville Formation was elevated to group status (Badgley, 1952) and 
correlations were extended into the subsurface (Badgley, 1952; Williams, 1963). A 
glauconitic sandstone found at the base of the Clearwater Formation was named the 
Wabiskaw Member of the Clearwater Formation (Badgley, 1952). 

The Mannville Group was subdivided into two units (Glaister, 1957 and 1959; 
Rudkin, 1964). Glaister (1957 and 1959) included the Wabiskaw Member with the upper 
unit, whereas Rudkin (1964) included it with the lower unit. 

The Blairmore Group of western Alberta was correlated with the Mannville Group 
on the basis of lithology, flora and fauna (Mellon and Wall, 1963). Correlations between 
the Athabasca and Lloydminster areas were made (Ruckin, 1964) and further stratigraphic 
work was undertaken by Mellon (1967) who commented on the difficulties in establishing 
the lower stratigraphic boundary of the Clearwater Formation on the basis of faunal 
changes. 

In the first geologic study of the Cold Lake area, correlations were made between 
the Coid Lake and Athabasca deposits (Webber, 1967). Correlations between all oil sand 


and heavy oil deposits were then made (Vigrass, 1968}. In the Lloydminster area, tongues 
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of marine sedimentary rock in the lower Upper Mannville Group were correlated to the 
Ciearwater Formation. In the Cold Lake area, the Mannville Group was divided into informal 
'A’,'B','C', and 'D' units (Vigrass, 1968) which were correlated approximately with the 
Upper Grand Rapids, Lower Grand Rapids, Clearwater, and McMurray Formations in the 
Athabasca area. The ’C’ unit top was placed slightly lower than the Clearwater Formation 
top in Athabasca. 

On the basis of three cores in the northeastern part of the Cold Lake area, 
Vigrass's (1968) informal 'A’ and 'B’ units were correlated with the Grand Rapids 
Formation, the 'C’ unit with the Clearwater Formation and the 'D’ unit with the McMurray 
Formation (Clack, 1968). 

Further refinements of correlation between all oil sand and heavy oil deposits were 
attempted (Kramers, 1974; Kendall, 1977; Mossop et al, 1981). These latter correlations 
underline the continuing problems with delimiting the upper and lower boundaries of the 
Clearwater Formation in the Cold Lake area, which will be discussed further below. 

The base of the Mannville Group is marked by an unconformity, while the Mannville 
Group-Colorado Group contact is considered unconformable (Mossop et al, 1981) or a 
diastem (Mossop, personal communication). 

The contact between the Clearwater Formation and the overlying Grand Rapids 
Formation has been described as a gradational interfingering facies boundary (Williams, 
1963) caused by time transgressive (Mellon and Wall, 1956-and 1963) or diachronous 
(Mellon, 1967; Mossop et al, 1981) deposition. 

The basal Clearwater contact with the underlying McMurray Formation is 
considered inter fingering (Mellon and Wall, 1956; Carrigy, 1963) or gradational (Mellon, 
1967). In northeast Alberta the change in the sandstone composition at the top of the 
McMurray Formation is still evident even though the sediments are faunally and 
lithologically gradational (Mellon, 1967). 

Petrographic criteria for differentiating between the McMurray and Clearwater 
Formations in the absence of fossils have been recognized (Carrigy, 1963). The 
Clearwater Formation can be recognized by the presence of glauconite pellets, euhedral 


biotite flakes, montmorillonite, and euhedral plagiociase grains. 
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B. Correlations derived from the Present Study 

Correlating the Clearwater Formation in the Cold Lake area is complicated by the 
fact that this area Is transitional between the northeastern Alberta oil sands deposits and 
the Lloydminster heavy oil deposit. Sedimentation Is significantly different in these areas. 
To the north, the Clearwater Formation is interpreted as an open marine shale, whereas in 
the south, its equivalent is interpreted as a complex mixture of marine shales, shelf and 
nearshore deposits, with some continenta! deposits. 

Previous attempts to correlate the Clearwater Formation in the Cold Lake area have 
adopted the stratigraphic nomenciature of the northeastern deposits (Figure 4). One 
notable exclusion is that the Wabiskaw Member is not included in the Cold Lake deposit 
since it cannot be distinguished here. 

The main correlation problem in the Clearwater Formation is the placement of the 
basal and upper shales. Previous studies (Kramers, 1974) placed the upper shale in the 
Lower Grand Rapids Formation and the basal shale with the McMurray Formation. More 
recently (Mossop, 1984), the upper shale has been included within the Clearwater 
Formation. 

To be stratigraphically and sedimentologically consistent, both the upper and basal 
shales should be included within the Clearwater Formation. In both the Athabasca and 
Wabasca oil sands areas, the McMurray / Clearwater and Clearwater /Lower Grand Rapids 
correlations are based on, in part, a lithologic change from the shale of the Clearwater 
Formation to the sands of the other two formations. In addition, there is no basis to 
exclude these shales from the Clearwater Formation in Cold Lake with respect to 
correlation with the Lloydminster area. The only micropaleontcologic sampie taken from the 
basal shale in the present study was poorly preserved, but was assigned on the basis of its 
dinoflagellate assemblage to the early Albian Age. Taking into account Mellon's (1967) 
comments on the difficulty of defining the base of the Clearwater Formation on faunal 
changes, the micropaieontologic results of the present study tend to confirm the 
placement of the basal shale within the Clearwater Formation. 

One basis for piacing the basa! shale within the McMurray Formation is that its 
contact with the Clearwater Formation is erosional. However, an erosional surface is often 


associated with the deposition of coarser-grained sediment. In addition, the nature of the 
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FIGURE 4. Stratigraphic correlation chart (modified from 
Kramers, 1975) 
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contact between the shale and the McMurray Formation is not Known since no core was 
recovered from this depth. If the basal contact of the shale is proven to be conformable 
with the McMurray Formation the proposed correlation could be incorrect. 

Further core and micropaleontologic work would help resolve this problem with 


stratigraphic correlations. 
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Ill. Sedimentology 


A. Previous Work 

The sedimentology of the oil sands and heavy oil deposits is highly compiex and a 
knowledge of facies is important to oil exploration (Mossop et al, 1981). 

Rudkin (1964) wrote extensively on the stratigraphy and depositional history of 
both the Athabasca and Lloydminster areas. !In his schema, the Lower Mannville Group 
included the McMurray Formation and the Wabiskaw Member of the Clearwater Formation 
and consisted of nonmarine basal fill deposits. The Upper Mannville Group included the 
Clearwater and Grand Rapids Formations and consisted of time transgressive marine units. 
The Clearwater Formation was considered to be a deltaic deposit, modified by periodic 
marine advances of a boreal sea whose axis trended southeast and lay close to the present 
Shield edge in northeastern Alberta. Marine conditions prevailed in the north (Athabasca 
area), mixed marine/non-marine conditions in the central area (Cold Lake and Lloydminster 
areas), and non-marine conditions in the south. 

The first description of the Clearwater Formation was made at the Athabasca 
outcrop area. The Clearwater Formation was described as interbedded sands and shales 
(McConnell, 1891) or as "soft grey shales, black shales, and grey and green sandstone 
with some hard concretionary layers” (McLearn, 1917). The first paleogeographic map 
(McLearn, 1932) depicted the Clearwater Formation as marine in origin. 

After further study, the Clearwater Formation was described as glauconitic 
sandstone basally, overlain by dark grey shales grading up into sandstone of marine origin 
(Mellon, 1967). Gallup (1974) studied the geologic history of McMurray-Clearwater 
deposition in the Athabasca oil sand area and characterized the environment of deposition 
as lagoonal and deltaic. 

A study of the chemical composition of the Mannville Group shales in central 
Alberta (Campbell and Williams, 1965) showed that Clearwater depositional environments 
were generally marine in the lower part and brackish in the upper part. 

A lithological description of the Mannville Formation in the Vermilion—Wainwright 
area was based upon the presence or absence of dark minerals, the distribution of marine 


fossils and the occurrence of coal seams (Nauss, 1945). The strata which are 
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approximately equivalent to the Clearwater Formation were described as: 88 feet of dark 
shales interbedded with silts and "salt and pepper” sandstone of marine origin; overlain by 
10 feet of quartz beach sandstone; overlain by 78 feet of deltaic "salt and pepper” 
sandstone, grey shale and siltstone with a few plant remains; overlain by 29 feet of beach 
reworked quartz sands, shale and siltstone. Overall, a nearshore to upper shoreface 
environment was indicated for the Vermilion—Wainwright area. 

Clearwater equivalent strata in the Lloydminster area (Wickenden, 1948) were 
described as composed of well-rounded fine-grained quartz with glauconite grains of 
probable delta front origin. Vigrass (1977) conducted a detailed study of the Lloydminster 
area and characterized the deposit as lithologically similar in all its nine members. The 
Clearwater Formation approximate equivalent Cummings to Sparky Members were 
considered marine in origin. Channel deposits in Clearwater equivalent Lloydminster strata 
(MacCallum, 1981) and nearshore marine deposits (Tilley and Last, 1980) were recognized. 

The 205 m (600 ft) of Mannville Group sediment in Coid Lake was described 
(Webber, 1967) as interbedded sand, silt and shale with minor stringers of coal. Ina 
detailed sedimentologic review of the Cold Lake area (Minken, 1974), the Clearwater 
Formation was characterized as a “salt and pepper” sandstone with glauconite, grading to 
shale southward. 

Cold Lake area sedimentology has been generally characterized as non-marine in 
the McMurray Formation, marine in the Clearwater Formation, and non-marine in the Grand 
Rapids Formation (Outtrim and Evans, 1977). 

The Clearwater Formation in the Cold Lake area has been generally described as a 
laterally extensive nearshore marine sand deposit with associated marine shales (Mossop 
etal, 1981). 

A general geologic report on the Clearwater Formation in the Cold Lake area 
(Carrigy and Kramers, 1974) revealed that marine conditions prevailed and that the Cold 
Lake deposit could be placed close to the southern limits of a boreal sea. Consequently, 
nearshore, transitional, and deltaic environments predominated. The Wabiskaw Member of 
the Clearwater Formation was described as shallow sea reworked McMurray Formation 


sands (Kramers, 1975). 
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The depositional environment of the Clearwater Formation at Cold Lake has been 
variously Characterized as nearshore deltaic (Jardine, 1974; Esso Resources, 1978), 
nearshore tidal (Wenneker et al, 1979), or an Arctic boreal sea transgressive marine bar 
and interbar sand wedge (Minken, 1974). 

In a site-specific study, Harrison et al (1981) examined over sixty wells in a three 
township area in the northeastern portion of the Cold Lake oil sands deposit. They 
concluded that the Clearwater Formation consisted of delta fringe silts and shales, lower 
delta front sands and shales, and upper delta front stream mouth bar and beach sands in 


this area. 


Provenance and Tectonic Framework 

Two source areas for Lower Cretaceous sediment have been proposed (Glaister, 
1959; Mellon, 1967; Carrigy and Kramers, 1974; Jardine, 1974; Keeler, 1980;, Harrison 
et al, 1981). The major source area was in the Cassiar-Omineca and Nelson uplift areas of 
south-central British Columbia, while minor contributions of sediment came from positive 
areas of the emergent Shield of Saskatchewan and Manitoba. The latter had a local 
significance in eastern Alberta, including the Cold Lake area. 

The hypothesis of dual provenance was reinforced by a geochemical study of 
Lower Cretaceous strata (Cameron, 1965) which found that the sediments could be 
divided into low, medium and high soda content groups which reflect different source 
areas. In the Athabasca, Cold Lake, and Lloydminster areas, tongues of low soda level 
sediment showed that some detritus came from the east, in addition to the major sediment 
source from the west. The Pre-Cambrian Shield to the northeast of the Lloydminster heavy 
oil deposit was cited as the source area for sediments deposited here by a southward 
transgressing Clearwater Sea (Putnam, 1979). However, in arecent paper, Putnam and 
Pedskalny (1983) contend that the source rocks for the Clearwater Formation at Cold Lake 
are only to the west and that there is no northeasterly Shield component. 

Relief on the Paleozoic erosion surface, which controlled Mannville Group 
deposition, was reported to be from 94 m (300 ft) to 120 m (400 ft) and locally up to 
150 m (500 ft) (Williams, 1963). An isopach and a paleogeologic map of the Lower 


Cretaceous strata in central Alberta north of the present area of study shows that rivers 
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flowed in a northwesterly direction in large channels etched into the pre-Cretaceous 
surface. The deposition of Mannville sediment was also controlled by a series of 
pre-Cretaceous islands trending in a northwest-southeast direction (Jardine, 1974). 
Mannville Group thickness was stated to be between 55U and 880 feet depending upon 
the sub-Cretaceous erosional relief (Clack, 1968). Two erosional pre-Cretaceous fluvial 
channels with northeriy drainage were recognized here. 

Williams and Stelck (1975) provided an overview of the depositional history and 
paieogeography throughout Mannville time in Alberta. A study of basement controls on 
Lower Cretaceous deposition (Stelck, 1975) showed that the Peace River and Sweetgrass 
arches controlled Lower Cretaceous sedimentation in Alberta. 

In the Cold Lake area, salt solution in late Mannville time, as well as subsequent to 
Mannville deposition, resulted in a discontinuity of facies (Energy Resources Conservation 
Board, 1973). Carbonate ridges and channels in the underlying Devonian strata also 
controlled sedimentation. 

Beaumont (1981) cites subsidence, caused by a combination of loading of thrust 
plates and erosion of these plates, as the major control on sea level and deposition in the 


Cretaceous Period. 


B. Introduction to Sedimentology 

Clearwater Formation is characterized by finely interbedded lithofacies. 
Consequently, it is most practical to describe the deposit with respect to facies 
associations. A facies association is a group of lithofacies which tend to occur together 
and are considered to be genetically or environmentally related (Reading, 1978). 

It is beyond the scope of the present study to perform a strict facies analysis of 
the cores. The main purpose of the study is to relate general sedimentological features 


and the geophysical well log response. 


Definition of Facies Associations 
The cores were studied bed-by-bed and sedimentary types were defined (See 


Table 1 for the range of features observed in core). 
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Table 1. 


Bed-by-bed 
Total Counts 


rae 


220 


151 


307 


140 


gay 


107 


39. 


105 


275 


426 


Totaly 267/ 


SEDIMENTARY TYPES 


Lithology 


Silty Shale 


Silty Shale 
Argillaceous Silt 


Carbonaceous Sand, 
Silt and Shale 
Coal 


Silty Shale, Shale, 
Argillaceous Silt 


Shale 


Sulait 
Silt 
Sand 


Calcareous Cement 


Sand 
Calcareous Cement 


Sand 
Calcareous Cement 


Sand 
Calcareous Cement 


Sand 
Sai kt 


All Lithologies 
All Lithologies 


Structures 


horizontal laminae; stringers; 
structureless 


laminated (horizontal, curved, 
inclined, undulating); load 
casts; wedges and lenses 
laminated (horizontal, curved, 
inclined, undulating); lenses; 
minor bioturbation 

heavily bioturbated 

laminated (horizontal, curved, 
inclined, undulating): load 
casts; stringers 


laminated (horizontal); 
stringers 


laminated (curved, inclined, 
undulating, contorted); lenses 


laminated (horizontal); 
cone-in-cone carbonate 


laminated (inclined); 
stringers 


laminated (undulating and 
curved); lenses and wedges 


apparently structureless 
heavily bioturbated 
remnant structures 
scoured surface 


silty shale rip-up clasts 
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The dominant lithologies in the Clearwater Formation are unconsolidated sand, 
unconsolidated silts to siltstone, silty shale, shale and calcareous cemented material. 
Silt-sized sediment is commonly oil saturated and unconsolidated and will be termed a silt. 
In the case of stratigraphic unit F, some silt-sized sediment is indurated and will be termed 
a siltstone. 

The major criteria for defining facies associations were grain size and the amount 
of bioturbation present. Heavy oil saturation in the sands and silts often obscures primary 
sedimentary structures and bioturbation often destroys these structures. 

The geophysical well log response was also used to help define the stratigraphic 
units and facies associations. It was the sole criterion used in situations where no core 


was available. The geophysical well log response is given in Appendix I. 


C. Facies Associations 

The Clearwater Formation is made up of six facies associations: heavily 
bioturbated silty shale; heavily bioturbated silty shale to silt and sand; moderately to heavily 
bioturbated sand with silt and silty shale; sand, silt and silty shale; graded sand, silt and 
shale with Bouma sequences; and sand and carbonaceous sand. 

The facies associations have been related to the stratigraphic units in the structural 


cross-sections A—B and A—C (Figure 5 to Figure 7). 


Facies Association 1. Heavily Bioturbated Silty Shale 

Facies Association 1 consists of a heavily bioturbated fissile silty shale (Plate 3-1). 
Some remnant horizontal laminae were preserved, although most primary sedimentary 
structures are destroyed by bioturbation. The silty shale is bioturbated by organisms, 


including those which produce Planolites nicholson (18793). 


Facies Association 2. Heavily Bioturbated Silty Shale to Silt and Sand 

Facies Association 2 consists of a grey shale with oil stained argillaceous silt and 
very fine-grained sand interbeds (Plates 1-1 and 1-2). The facies association becomes 
more shalely to the northwest. Silty shale thicknesses range from individual laminae to 


beds 0.45 m (1.5 ft) thick. Silt and sand range from individual laminae to beds up to 20 cm 
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silty shale interlaminations 


coal interlaminations 


calcareous cement 


| argillaceous silt interlaminations 
i—-— 
vee 


«K 


cone-in-cone_ structure 


parallel laminations 

inclined laminations 

ripple cross-laminations 
bi-directional cross-laminations 
erosional surface 


syndepositional fault 


Micropaleontology 


Sample location 


Brackish 


Nearshore Marine 


Open Marine 


horizontal burrows 


vertical burrow 


rip-up clast 


carbonaceous debris 


macrofossil 


correlation based on core and geophysical well logs 


correlation based partialiy or entirely on geophysical well logs 


(core missing) 


FIGURE 5. Legend for structural cross-sections. 
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(8 in) thick. The thickness of sand and silt interbeds increases upwards within stratigraphic 
units J to |. 

The beds typically consist of horizontal and undulating laminae with rare inclined 
cross-laminae dipping at a low angie of 4°. Some bi-directional low-angle cross-laminae 
dipping at 8° were observed in the basal portion of the facies association in well 
10-7-66-5W4. Small silty shale clasts of iess than 1 cm diameter "float" in some of the 
oil saturated silt and sand beds. Synsedimentary normal faults were observed in well 
10-17-66-5W4. They are not coring artifacts because they occur within the core and do 
not extend to the core edge. 

The shales and silty shales are highiy bioturbated by traces of Planolites nicholson 
(1873). The location of bioturbation is given in Figures 5 to 7. The degree of bioturbation 


decreases upwards within stratigraphic units J to |. 


Facies Association 3. Moderately to Heavily Bioturbated Sand with Silt and Silty 
Shale 

Facies Association 3 consists of interbedded to interlaminated very fine-grained 
sand to silt and silty shale (Plates 1-3, 2—1 and 3-1). Sand units range from individual 
laminae to beds 0.9 m (3 ft) thick in well 6-3-66-5W4. Silty shale is found as individual 
laminae to beds 25 cm (10 in) thick. The proportion of silty shale increases towards the 
northwest. Carbonaceous debris is found in a 10 cm (4 in) thick interbed of laminated silt in 
well 6-3-66-5W4. In stratigraphic unit G, carbonaceous debris is found in silty shale 
beds up to 4 cm (1.5 in) thick. 

Structures are obscured in the heavily oil saturated sands and in stratigraphic 
subunit D, where bioturbation has destroyed primary sedimentary structures. Where 
observed, structures in the sands and silts are dominately horizontal laminae. Undulating 
and ripple cross- laminae are also commonly found. Some sequences (up to 25 cm thick in 
stratigraphic unit H in well 6-3-66-5W4) of interbedded and interlaminated sand and silty 
shale are inclined up to 15° and are probably accretion surfaces. Bi-directional low-angle 
inclined laminae were noted in stratigraphic unit G in well 2—19-66-5W4 and 
bi-directional ripple cross-laminae were noted in stratigraphic unit G in well 


6-5-66-5W4. 
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The silty shale exhibits horizontal and inclined laminae. The low-angie inclined 
laminae dip at angles between 3° and 10°. The silty shale laminae often exhibit scoured 
surfaces. 

Silty shale intraclasts (0.1 cm to 2 cm in diameter) are numerous in some parts of 
the core "floating’ within the oil saturated sands. The silty shale intraclasts are most likely 
derived from silty shales which have undergone scouring. 

Bioturbation of the sands, silts and silty shales is extensive in stratigraphic unit D, 
and the basal portion of stratigraphic unit H. The degree of bioturbation in the silty shales 
decreases upwards through stratigraphic units H and G. Bioturbation consists of 1 cm 
diameter traces of Planolites nicholson (1873) and simple vertical traces of Skolithos 
linearis haldeman (1840). The latter are found in association with Planolites nicholson 
(1873) at equivalent depths of stratigraphic unit H in wells 6—5-66-5W4, 
10-—7-66-5W4 and 10-—17-66-5W4. This trace varies between 1 cm and 2 cm in 
diameter and cuts across sedimentary types, attaining a maximum length of 22 cm in well 


10-17-66-5W4. 


Facies Association 4. Sand, Silt and Silty Shale 

Facies Association 4 consists of very fine-grained, moderately poorly sorted to 
poorly sorted sand, interlaminated or interbedded with poorly sorted silty shale, 
argillaceous silt and silt (Plates 2-2 and 2-3). The proportion of silty shale laminae and 
beds increases to the northwest in stratigraphic subunit F,. 

Sand and silt are found primarily as horizontal laminae, but beds can attain 
thicknesses of up to 20 cm (6 in). The silty shale is also usually finely interlaminated with 
the coarser units, although silty shale beds up to 15 cm (6 in) thick occur. Minor portions 
of the sands and siltstones in stratigraphic subunit F, exhibit low-angle cross-laminae 
inclined up to 7°. To the northwest, some small (up to 0.5 cm diameter) silty shale 
intraclasts were noted in this facies association. A small synsedimentary normal microfault 
was noted in stratigraphic unit E in well 2-19-66-5W4. The fault does not extend to the 
core edge and the core does not exhibit features suggestive of breakage. 

Bioturbation is rare and consists of small (up to 0.5 cm diameter) isolated traces of 


Planolites nicholson (1873) in the silty shales. in the upper silty shale bed of stratigraphic 
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subunit F, in well 10-17-66-5W4, a shell specimen of Psilyma peterpondia and a trace 
of Muensteria von Sternberg (1883) were found. 

Calcareous cemented siltstones with associated cone-in-cone structures were 
observed in stratigraphic subunit F, in wells 2-19-66-5W4 and 10-17-66-5W4. The 
cone-in-cone structures suggest that high pressure levels existed during dewatering of the 


sediment. 


Facies Association 5. Graded Sand, Silt and Shale with Bouma Sequences 

Facies Association 5 consists of interbedded to interlaminated sand, silt and silty 
shale. 

The sand is found as individual laminae and as beds up to 0.6 m (2 ft) thick (Plate 
3-2). Silty shale is found as individual laminae and in beds of laminated silty shale up to 13 
cm (5 in) thick. 

Major structures include undulating, horizontal and low-angle inclined (up to 10°) 
laminae. Bi-directional low-angle cross-laminae were observed in all wells. The silty shale is 
commonly scoured. Silty shale intraclasts are commonly found. These clasts are 
subangular and average |.5 cm in diameter. Sequences in well 10—-17-66-5W4 between 
443.5 m (1454 ft) and 444.5 m (1457.5 ft) exhibit graded beds with Bouma Sequences 
(Plate 3-2). These deposits consist of a very fine-grained, heavily oil saturated, massive 
sand (Bouma Ta subdivision), grading upwards into a lightly oil stained, plane parallel 
laminated silt (Bouma Tb subdivision). Intraclasts are incorporated into either the sand or 
silt beds. The sequence is capped by ripple cross-laminated sands, draped by undulating 
silty shale laminae (Bouma Tc subdivision) and parallel silty shale laminae (Bouma Td 
subdivision). Some of the horizontally laminated silty shale represent interturbidite 
deposits (Bouma Te subdivision). The topmost laminae of shale is usually scoured, marking 
the commencement of the next turbidity current deposit. Rare Planolites nicholson (1873) 


traces were noted. 
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Facies Association 6. Sand and Carbonaceous Sand 

Facies Association 6 consists of very fine-grained, poorly sorted sand with rare 
silty shale laminae and silty shale stringers. Carbonaceous debris is concentrated in some 
silty shale laminae (Plates 2-2 and 3-3). 

Sand is found as individual laminae or as large beds of laminated sand up to 12 m 
(40 ft) thickness in stratigraphic subunit F,. Silty shale is found as individual laminae or as 
beds of laminae up to 5 cm (2 in) thick in stratigraphic unit B. Carbonaceous debris and a 
coal lamina 2 cm (0.8 in) thick were observed in the basal portion of stratigraphic unit B in 
well 10-17-66-5W4. 

Sedimentary structures within this facies association are predominately horizontal 
laminae, inclined low-angle cross-laminae (up to 12°) and ripple cross-laminae. The ripple 
cross-laminae are commonly draped by undulating silty shale laminae. Bi-directional 
low-angle cross-laminae were observed in stratigraphic subunit F, in well 
10-17-66-5W4. Scoured surfaces were seen in stratigraphic subunit F, in most wells. 
Usually the scouring is associated with individual or groups of silty shale intraclasts up to 2 
cm in diameter. A small synsedimentary normal microfault is associated with an erosional 
surface in stratigraphic subunit F; in well 6-3-66-5W4. The fault does not extend to the 
core edge and is not a coring artifact. Rare individual traces of Planolites nicholson (1873) 
occur in Facies Association 4. 

Small portions (up to 15 cm (6 in) thick) of stratigraphic subunits F, and F; are 
altered by calcareous cementation, but some remnant ripple cross-laminae and horizontal 


laminae are preserved. 


D. Stratigraphic Units 

The Clearwater Formation was divided into ten informal stratigraphic units on the 
basis of lithology and geophysical well log response (Figures 5 to 7 and Appendix |). The 
stratigraphic units were defined to describe groups of facies associations which could be 
combined into higher order assemblages. The informal stratigraphic units were also used 
in the mapping and delineation of large scale vertical sequences to define regional 
sedimentation patterns. This required an ordering of the system at a level higher than that 


of facies associations. 
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Stratigraphic unit J (Facies Association 2) overlies the McMurray Formation 
sediments. The nature of the contact between the McMurray and Clearwater Formations in 
this area is unknown since no core was available. The thickness of stratigraphic unit J 
ranges from 2.5 m (8 ft) to 5.5 m (18 ft). Core from this unit is available from only three 
of five wells, with a maximum length of 1.5 m (5 ft) cored in well 6-5-66-5W4. 
Stratigraphic unit H exhibits a gamma ray response of 128 API units and a deep induction 
log response of 2 ohms m?/m. 

Stratigraphic unit | (Facies Association 2) overlies stratigraphic unit J with either a 
sharp or erosional contact. In wells 6-5-66-5W4 and 10-7-66-5W4 the contact is 
sharp, but cannot be conclusively identified as erosional. However, well 10—17-66-5W4 
has an erosional contact. The thickness of this stratigraphic unit is fairly uniform, ranging 
from 2.3 m (7.5 ft) in the southeast to 3.8 m (12.5 ft) in the northwest. Complete core 
from this stratigraphic unit was available in three out of five wells (6-5—66—-5W4, 
10-7-66-5W4 and 10-17-66—-5W4). Stratigraphic unit | was recognized from the 
geophysical well logs by an abrupt shift in the gamma ray log from 128 API units to 95 API 
units at the top of stratigraphic unit |. The deep induction log value for this stratigraphic 
unit is 3.5 ohms m?/m. 

Stratigraphic unit | is overlain by stratigraphic unit H (Facies Association 3). The 
nature of the contact is usually gradational, but in the case of well 10-—17-66-5W4 it is 
sharp. The differentiation of stratigraphic units | and H is based upon a change in the 
lithology and in the degree of bioturbation. In the case of a gradational contact, the degree 
of bioturbation in the silty shales is similar between stratigraphic units | and H and the 
contact must be picked from geophysical well logs. Stratigraphic unit H has a variable 
thickness, ranging from 3.7 m (12.5 ft) in the northwest and 5.2 m(17.5 ft) in the 
southeast to 6 m (20 ft) at well 1O—7-66-—5W4. Complete core was available from three 
out of five wells. In addition, well 6-3-66—-5W4 is missing only the basal 0.8 m (2.5 ft) 
and well 2—-19-66-5W4 has only the upper 0.6 m (2 ft) cored. The base of stratigraphic 
unit H was picked at the 95 API unit level on the gamma ray log. The deep induction log 
response increases uniformly upward from 3.5 ohms m?/m at the base to 8 ohms m?/m 


at the top of this stratigraphic unit. 
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Stratigraphic unit G (Facies Association 3) gradationally overlies stratigraphic unit H. 
The contact must be picked with the aid of geophysical well logs. In well 6-5-66-5W4, 
stratigraphic unit G is less well developed, but the contact with stratigraphic unit H is more 
apparent. Thickness of stratigraphic unit G sediments is variable, ranging from 0.75 m (2.5 
ft) in well 6-5-66-5W4 to 3 m(10 ft) in well 6-3-66-5W4 to 4.5 m (15 ft) in the other 
three wells. Core recovery is good with only small parts of core missing. The contact 
between stratigraphic units G and H is marked by a deflection of the gamma ray log value 
from 90 API units to 83 API units. The deep induction log displays a sharp increase from 8 
ohms m?/m to 12 ohms m?/m. 

Stratigraphic unit F sharply overlies stratigraphic unit G. it is the thickest unit in the 
Clearwater Formation, averaging 13.75 m (45 ft) in the five wells studied. Core recovery 
is fair with minor portions missing. Stratigraphic unit F has been divided into three 
subunits, primarily on the basis of lithology. Stratigraphic subunit F; (Facies Association 6), 
the basal portion of stratigraphic unit F, ranges in thickness from 4.2 m (14 ft) in the 
northwest to 13 m (42 ft) in the southeast. The contact between stratigraphic subunit F, 
and the underlying stratigraphic unit G is marked by a deflection of the gamma ray from 87 
API units at the top of stratigraphic unit G to 75 API units at the base of stratigraphic 
subunit F,. The deep induction log registers an increase from 11 ohms m?/m at the top of 
stratigraphic unit G to 13 ohms m?/m in stratigraphic subunit F,. Stratigrahpic subunit F, 
(Facies Association 4) thickens to the northwest. It ranges in thickness from 6.4 m (21 ft) 
in well 2—-19-66-—5W4 to 2.4 m (8 ft) in well 6-3-66-5w4. Core recovery was good in 
the sands, but only fair in the shales. Some shale core was destroyed during coring and 
correlations had to be established through an examination of the geophysical well logs. 
Stratigraphic subunit F, is distinguished from stratigraphic subunit F; by an increase in the 
gamma 75 API units to 90 API units. The deep induction log decreases abruptly from 13 
ohms m?/m to 6 ohms m?/m. The thickness of stratigraphic subunit F, (Facies Association 
6) is highly variable. The thickest section is 3 m (10 ft) and is found in well 
10-7-66-5W4. The thinnest section is 0.75 m (2.5 ft) and is found in well 
6-3-66-5W4. Core recovery is poor. Stratigraphic subunit F, is distinguished from 
stratigraphic subunit F, by a deflection of the gamma ray log value from 97 API units at the 


top of stratigraphic subunit F, to 77 API units near the base of stratigraphic subunit F,. The 
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deep induction log increases rapidly to 30 ohms m?/m in stratigraphic subunit F, from 9 
ohms m?/m at the top of stratigraphic subunit F,. 

Stratigraphic unit E (Facies Association 4) sharply overlies stratigraphic unit F. It has 
a variable thickness, ranging from 3.4 m (11 ft) in well 6-3-66-5W4 to 0.9 m (3 ft) in 
well 10—7-66-5W4. The base of stratigraphic unit E is marked by a gamma ray value of 
98 API units and a rapid decrease in the deep induction log to 11 ohms m?/m. 

Stratigraphic unit D sharply overlies stratigraphic unit E. This stratigraphic unit 
thickens to the northwest. Thicknesses range from 2.1m (7 ft) in well 6-3-66-5W4 to 
4.5m (15 ft) in well 2—-19-66-—5W4. Core recovery was poor. Wells had to be correlated 
using geophysical well logs. Stratigraphic unit D is made up of three subunits. !t consists 
of avery fine-grained, poorly sorted sand unit (stratigraphic subunit D,) enclosed by fissile 
silty shales (stratigraphic subunits D, and D,). Stratigraphic subunit D, (Facies Association 1) 
ranges in thickness from 0.5 m (1.7 ft) to 1.75 m (5.8 ft) in well 10—17-66-5W4. It is 
distinguished from stratigraphic unit E by a deflection of the gamma ray log to 108 API 
units and a decrease in the deep induction log response from 11 ohms m?/m to 8.5 ohms 
m?/m. Stratigraphic subunit D, (Facies Association 3) ranges in thickness from 0.5 m(1.7 
ft) in the southeast to 2 m (6 ft) in the northwest. The gamma ray log value shifts slightly 
from 108 API units to 102 API units and the deep induction log exhibits a sharp increase 
from 8.5 ohms m?/m to 13.5 ohms m?/m. Stratigraphic subunit D, (Facies Association 1) 
ranges in thickness from 1 m (3.3 ft) to 2.75 m (9.1 ft) in well 10—17-66—-5W4. it is 
distinguished from stratigraphic subunit D, by an increase in the gamma ray log value from 
102 API units to 108 API units and by a deep induction log response decreasing from 11 
ohms m?/m to 7 ohms m?/m. 

Stratigraphic unit C (Facies Association 5) sharply overlies stratigraphic unit D. The 
thickness of this stratigraphic unit is highly variable, ranging from 2.7 m (9 ft) in well 
10-17-66-5W4 to 6.2 m (21 ft) in well 6-5-66-5W4. In general, this stratigraphic unit 
thickens to the southeast. Core recovery is good except in well 2—- 19-66-5W4, where 
only 2.1 m(7 ft) out of 3.7 m (12 ft) was cored. Stratigraphic unit C was recognized on 
the gamma ray log by a sharp decrease from 108 API units to 87 API units and on the deep 


induction log by an increase from 9 ohms m?/m to 13 ohms m?/m. 
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Stratigraphic unit B (Facies Association 6) sharply overlies stratigraphic unit C. It is 
fairly uniform in thickness, averaging 1.5 m (5 ft). Core recovery is poor. Well 
10-—17-66-5W4 was the only well from which complete core was recovered. 
Stratigraphic unit B is characterized by a gamma ray response of 82 API units and a rapid 
increase in the deep induction log from 12 ohms m?/m to 25 ohms m?/m. 

There was no core recovery from stratigraphic unit A. It is recognized by an 
increase in the gamma ray response from 101 AP! units to 117 API units and a marked 
decrease in the deep induction log from 12 ohms m?/m to 4 ohms m?/m. These 
responses suggest that stratigraphic unit A is a shale. The thickness of this stratigraphic 


unit ranges from 4.5 m (15 ft) in the northwest to 6 m (20 ft) in the southeast. 


E. Micropaleontology 

A total of fourteen samples were submitted to Dr. Singh (Alberta Geological! 
Survey) for analysis of dinoflagellate and megaspore assemblages. The location of the 
samples is given in the structural cross-sections (Figure 6 and Figure 7). 

Eighteen dinoflagellate species were identified. Four samples were poorly 
preserved and must be interpreted with caution. The moderately preserved samples 
contained up to 20 specimens. There were no well preserved samples. 

Dinoflagellate populations are useful in stratigraphic studies of the Cretaceous 
since they give highly accurate stratigraphic results (Sarjeant, 1974). All samples submitted 
were assigned to the early Albian Age. 

The use of dinoflaggelate assemblages in paleoenvironmental reconstruction is less 
conclusive since microfossils are easily transported in shallow marine environments 
(Heckel, 1972). Paleoenvironmental interpretation of dinoflagellate assemblages is based 
on population diversity (C. Singh, personal communication), species morphology, and 
absolute frequency (Sarjeant, 1974). The presence of few species suggests low salinity 
brackish conditions. Thin-walled, less ornamented cysts are found in shallow, more turbid 
environments, while more delicate specimens denote deeper, less turbid conditions. An 
increase in cyst density indicates open marine conditions. But as the proportion of 
land-derived material declines, total numbers drop. Thus fewer specimens can be 


indicative of either nearshore or deep-water conditions. The seaward extension of river 
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deltas can create conditions which change the dinoflagellate assemblage. The large 
amounts of detritus in the water column would decrease the light supply (Sarjeant, 1974) 
and the assemblage would reflect a brackish environment. 

An examination of samples for the five wells in this study confirms the complexity 
of sedimentation in a shallow epeiric sea. 

Stratigraphic unit J was sampled in well 10-17-66-5W4. Recovery and 
preservation were poor, but the presence of one specimen of Mudergonia assymetrica 
(Brideaux, 1977) and another specimen of Mudergonia sp. suggests brackish conditions. 
The specimen in stratigraphic unit J is from the basal shale. The existence of brackish 
conditions at the time this shale was deposited suggests that fully open marine conditions 
did not exist in this area of the shallow sea. Further investigation is recommended to 
confirm this interpretation, especially because the interpretation is based on only two 
specimens. The poor preservation of these specimens suggests that they may have been 
transported rather than part of the “living community”. 

The number of samples examined is too low for regional! interpretations to be 
made, but a general trend toward open marine and nearshore marine envircnments is 
reinforced by paleoenvironmental interpretations of stratigraphic units H and F. The 
assemblage in stratigraphic unit F (well 2-19-66-5W4) is consistent with a more saline 
environment. However, the presence of nine specimens of megaspores, comprising 
Bacutriletes, Erlansonisporites, and Minerisporites suggests proximity to a shoreline which 
was not recognized in core. 

Examination of specimens in Facies D suggest that brackish conditions prevailed 
towards the northwest and open marine conditions prevailed towards the southeast 
portion of the study area. However, in an area of subsidence such as this, low salinity 
brackish conditions can prevail just before open marine conditions are established (Singh, 
personal communication). These species then die out when open marine condistions are 
fully established and are replaced by a suite of open marine assemblages. Stratigraphic unit 
D in wells 2—-19-66—-5W4 and 10-17-66-5W4 reflects this pattern. Samples taken at 
the base of stratigraphic unit D are brackish. Preservation was poor in this sample. 
However, 9 dinoflagellate specimens were identified along with thirteen megaspore 


specimens, comprising Erlanisporites and Bacutriletes. Further upsection, at 448.9 m 
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(1459.5 ft) in well 10-17-66-5W4, open marine assembiages are found. Samples taken 
from the top of stratigraphic unit D and from the overlying stratigraphic unit C contain 
poorly preserved brackish and nearshore suites, suggesting a return to more less saline 
conditions. Once again, no actual evidence for emergence was found. This would suggest 
that the shoreline shifted rapidly and was subjected to reworking in this shallow sea 
environment. 

Micropaleontologic data reflect deposition in a "shelfal” to nearshore environment. 
More detailed sampling, along with sedimentologic study is required to confirm 


interpretations in a complex environment such as this. 


F. Sedimentary Textures 

The textures of the sediments in the Clearwater Formation were analyzed using 
standard sieve analysis procedure. Results are given in Appendix Il. 

A plot of standard deviation (sorting) versus mean grain size (Figure 8) shows that 
the sands are uniformly very fine-grained (between 3 and 4 phi) and predominantly poorly 
sorted (between 1 and 2 phi standard deviation according to Folk's (1968) classification) 
(Plates 6-1 and 6-2). An excellent direct relationship (r=0.94) is exhibited between the 
standard deviation and the mean grain size. The coarsest samples are the best sorted, 
whereas the finer-grained samples are the most poorly sorted. 

General features of sand, argillaceous sand and matrix-supported samples are 
shown in representative histograms. Relatively clean sands (Figure 9) are unimodal and 
show a pronounced peak at the 4 phi level. There are relatively few fines in these samples. 
Argillaceous oil sands (Figure 10) display the same unimodal peak as the oii sands samples, 
but there is amore pronounced skewness towards the fines (higher phi values). 
Matrix-supported samples exhibit bimodality (Figure 11). The major peak is found at the 
3.5 to 4 phi level which is similar to previously-mentioned sample types. In addition, a 
more diffuse peak is located between 9 phi and 10.5 phi. This relects what is seen in thin 
section as discrete grains “floating” in a fine matrix. 

Two coarsening-upward sequences can be discerned within the very fine-grained 
sands, one from stratigraphic unit | to the top of stratigraphic unit F and the other from 


stratigraphic unit E to stratigraphic unit B (Figure 8). The sands in this latter 
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FIGURE 8. Standard deviation (sorting) versus mean grain size, 
showing relationship of stratigraphic units to grain 
size parameters. (r = 0.94) 
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FIGURE 9. Representative grain size histogram for 
Sand samples. 
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coarsening-upward sequence are slightly coarser-grained, although they are still within the 


very fine-grained classification. 


G. |lsopach and Structural Mapping 

The geophysical well log reponse for each stratigraphic unit (Appendix |} was used 
to correlate facies regionally. Fifty-four wells (Figure 12) were examined in a ten township 
area (Townships 64 to 66, Ranges 4W4 to 6W4 and Township 64, Range 3W4). 
Structural maps were generated for the pre-Cretaceous unconformity and the tops of 
stratigraphic units F, G, H and A to discern the regional dip of these features. lsopach 
maps were generated for all stratigraphic units to delineate regional sedimentation trends 
and to examine the effects of the structure of the pre-Cretaceous unconformity on 


subsequent sedimentation. 


Structural Maps 

The structural maps of the pre-Cretaceous unconformity (Figure 13) agrees well 
with the structural map of the pre-Cretaceous unconformity constructed by Martin and 
Jamin (1963) and Kramers (1975). The pre-Cretaceous surface dips to the south-southeast 
at arate of about 5 metres per kilometre. A local high of the pre-Cretaceous 
unconformity trending northwest-southeast is found in the western part of Township 
66-5W4, the eastern part of Township 66—6W4 and the northeastern part of Township 
65-6W4. Another local high trending northeast-southwest is located in Township 
64-6W4. A local low is found in the southern part of Townships 65-5W4 and 64-5W4. 

The structural maps of the tops of stratigraphic units F, G, H and the top of the 
Clearwater Formation (stratigraphic unit A) depart from that of the pre-Cretaceous 
unconformity, but are similar to one another. This indicates tectonic movement prior to 
Clearwater Formation deposition and a uniform response to tectonic movement during the 
deposition of Clearwater sediments. 

The deposits of the Clearwater Formation dip from north to south at arate of 2.5 
metres per kilometre. In the central part of the area there is a structural high trending 
northwest-southeast. A structural low trending northeast-southwest is located in 


Township 64-—5W4. These features are found in all facies examined, but are best 
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FIGURE 12. Regional well control for structural 
and isopach maps. 
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FIGURE 13. Structural map of pre-Cretaceous unconformity. 
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So 
expressed in the structural map of the top of stratigraphic unit F (Figure 14). 


isopach Maps 


Stratigraphic Unit J 

The isopach of stratigraphic unit J (Figure 15) shows a blanket shale which 
attains a maximum thickness of 6 m (20 ft). A north-south trending thinning is 
located in the central part of the region where thicknesses are of the order of 0.9m 
(3 ft) to 2 m (7 ft). The distribution of thicknesses is very dissimilar from other 
stratigraphic units in the Clearwater Formation. This is further evidence for the 
placement of an erosional surface between stratigraphic units J and |. In 
stratigraphic unit J there is a small area of thickening is the northwest corner of 
Township 65—-5W4 situated over one of the local pre-Cretaceous highs. Here the 
stratigraphic unit reaches a thickness of 6 m (20 ft). This represents the most 
complete section of this stratigraphic unit in the area. It may have been protected 
from erosion by its location on the local pre-Cretaceous high, above the general 


level of erosion. 


Stratigraphic Unit | 

Stratigraphic unit | (Figure 16) is thickest in a line trending 
northeast-southwest where it is up to 7.2 m (24 ft) thick. The pre-Cretaceous high in 
Township 64-—6W4 has some control on the development of a thicker deposit here, 
but the other highs do not have an effect on sedimentation. The thinnest 
accumulations of stratigraphic unit | are between 1.5 m (5 ft) and 2.1 m (7 ft) and are 


found in the extreme northwest and southeast of the area. 


Stratigraphic Unit H 

Stratigraphic unit H (Figure 17) is almost the mirror image of stratigraphic unit 
|. The sediments of stratigraphic unit H are up to 10.5 m (34 ft) thick in in the 
northwestern and southeastern parts of the area and thinnest in an area trending 
northeast to southwest through the centre of the region. Here the stratigraphic unit 
is as thin as 1.8 m (6 ft). The local high in Township 64-—6W4 still exerts control on 


sedimentation and a thickening of the stratigraphic unit to 8.2 m (27 ft) is seen here. 
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FIGURE 14. Stratigraphic unit F structural map. 
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FIGURE 15. Stratigraphic unit J isopach map. 
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FIGURE 16. Stratigraphic unit I isopach map. 
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FIGURE 17. Stratigraphic unit H isopach map. 
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Stratigraphic Unit G 

Stratigraphic unit G (Figure 18) reaches a maximum thickness of 4.5 m (15 ft) 
over the pre-Cretaceous high in Township 64-6W4. It is thinnest in an east-west 
trending line through the central part of the area where it is 1.2 m (4 ft) to 2.1 m (7 


ft) thick. 


Stratigraphic Unit F 

The stratigraphic unit F isopach map (Figure 19) reveals a departure from 
previous sedimentation patterns. The thickest deposits are found ina ine trending 
southeast-northwest through the centre of the region. Here accumulations of 
between 15.5 m (51 ft) in the northwest and 26.5 m (87 ft) in the southeast are 
found. The thick accumulation of sand in the northwest part of Township 64-5W4 
and the southwest part of Township 65-5W4 represents the formation of large 
scale features. 

The thinnest accumulations are 9.5 m (32 ft) and are located in the northeast 
and southwest parts of the region. The pre-Cretaceous unconformity has little 
effect on the sedimentation pattern for stratigraphic unit F as a whole. 

The sand unit enclosed by silty shale (stratigraphic subunit F,} was also 
studied for regional patterns (Figure 20). The silty shale was identified on the basis 
of an increase in the gamma ray response as compared with the rest of the 
stratigraphic unit. The actual gamma ray API cutoff varied between wells. Only one 
silty shale unit could be identified in some wells indicating that this stratigraphic 
subunit is not regionally correlative. The total thickness of silty shale increases from 
0.6 m (2 ft) in the south-central region to a high of 9 m (30 ft) in the north and west. 

Thinnest accumulations of shale in stratigraphic subunit F, occur towards the 
southeast, while thicker accumulations occur towards the north and west, 
suggesting more open marine conditions prevailed in these areas. 

The sand/ shale ratio map (Figure 21) for this portion of stratigraphic unit F 
shows that the areas with thicker shales also have a thicker sand unit. Sand/ shale 
ratios range from O in areas where there was only one shale bed to 0.75 in the 
British Petroleum Pilot Plant area (Township 66-5W4). There are definite loci of 


greater sand content. Some of these areas correspond to highs on the 
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FIGURE 18. Stratigraphic unit G isopach map. 
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FIGURE 19. Stratigraphic unit F isopach map. 
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FIGURE 20. Stratigraphic subunit F, isopach map. 
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FIGURE 21. Stratigraphic subunit Fo sand/shale ratio map. 
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pre-Cretaceous unconformities, whereas others do not. The patches of higher 
sand/ shale ratio values in the southeast part of the region reflect major sites of 


sand deposition. 


Stratigraphic Unit E 

Stratigraphic unit E is very thin and was recognizable in only a few of the 
more recently logged wells. It has a very erratic distribution in the northwest part of 
the region and was not mapped. Its limited extent may be a result of the limitation of 
the geophysical well log or may reflect an actual limited distribution of this 
stratigraphic unit, suggesting localized processes for the deposition of this 


sediment. 


Stratigraphic Unit D 

The distribution of stratigraphic unit D is also erratic (Figure 22). It is thickest 
in a line trending northeast-southwest in the southwestern part of the region. Here it 
attains a maximum thickness of 7.2 m (24 ft). This stratigraphic unit is locally thick in 
the northwestern part of the region (Township 66—6W4) where it is 6.6 m (22 ft) 
thick. The highs on the pre-Cretaceous unconformity have some influence on the 
development of thicker accumulations of stratigraphic unit D in this region. This 
stratigraphic unit has a greater thickness over a high in the southeastern part of the 


region where accumulations reach 6.4 m (21 ft). 


Stratigraphic Unit C 

Stratigraphic unit C (Figure 23) continues the trends of stratigraphic unit D 
and attains thicknesses of 5.4 m (18 ft) in the southwest and 6.4 m (21 ft) in the 
northwest. The effect of the pre-Cretaceous unconformity in influencing thicker 
accumulations is more pronounced in the southwest. The thinnest accumulation is 
1.8 m (6 ft) found in the western part of the region trending northeast-southwest. 
Stratigraphic unit C is also thick in the south-central part of the region where it 


attains a thickness of 4.8 m (16 ft). 
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FIGURE 22. Stratigraphic unit D isopach map. 
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FIGURE 23. Stratigraphic unit C isopach map. 
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Stratigraphic Unit B 

The distribution of stratigraphic unit B (Figure 24) is similar to that of 
stratigraphic units C and D. Accumulations in the northwest and southwest are up to 
2.1m (7 ft) over the highs of the pre-Cretaceous unconformity and up to 2.4 m (8 
ft) in the east-central area. Thinnest accumulations are again found in a 


northeast-southwest trend where only 0.6 m (2 ft) of stratigraphic unit B is present. 


Stratigraphic Unit A 

Stratigraphic unit A (Figure 25) continues the distribution trend, although it is 
less variable in thickness. The thickness of this stratigraphic unit ranges from 7.7 m 
(26 ft) in the northwest and southwest to 6.6 m (22 ft) in the south-central portion 
of the region to 4.5 m (15 ft) in a northeast-southwest trend in the central part of 


the region. 


H. Facies Association Interpretation 


Introduction 

The mineralogy throughout the Clearwater Formation in this area is a feldspathic 
litharenite to litharenite (Folk, 1968) with considerable glauconite. Texturally, the 
sediments are immature to submature (Folk, 1968). The abundance of glauconite and 
absence of any subaerial features (i.e., dessication cracks, rooted zones, caliches, etc.) 
suggest that the environment of deposition was marine. 

The facies associations and the succession of stratigraphic units in the British 
Petroleum Pilot Plant area are the result of outer and inner shelf deposition in a shallow 
marine environment which received very fine-grained detritus from the southeast. 

The facies association interpretation and summary of stratigraphic units (Figure 26) 
is applicable only to the Clearwater Formation in the Cold Lake oil sands deposit and is a 
local summary rather than a model (Walker, 1979a). 

The study of shallow epeiric seas is complicated by limitations in terminology that 
can be applied to the environments of deposition. The gradient of the epeiric sea floor is 
low. No discernable s!ope break is evident in shallow epeiric seas and there is a gradual 


transition to deeper water sedimentation. 
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FIGURE 24. Stratigraphic unit B isopach map. 
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FIGURE 25. Stratigraphic unit A isopach map. 
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Top of Clearwater Formation 
A Shate, not cored 
B Carbonaceous sand, minor silty shale; 


laminated 
Carbonaceous sand with silty shale; 


Silty shale and sand; interbedded to 
interlaminated; heavy bioturbation 
Sand and silt; finely horizontally 


__ silty shale at top 

Sand; horizontally laminated, cross- 
_ and minor ripple cross-laminated 
Sand and silty shale; horizontally 
laminated, cross and minor ripple 
cross-laminated; minor bioturbation 
Sand with minor silty shale; horizon- 
tally laminated, cross- and minor 
ripple cross-laminated; minor 
bioturbation; erosional surfaces 

and rip-up clasts 

Carbonaceous sand with silty shale; 
interbedded, cross-laminated; minor 
bioturbation; rip-up clasts 

Sand with silty shale; horizontally 
interlaminated to interbedded; 
moderate bioturbation 

Sand and silty shale; horizontally 
interlaminated to interbedded; 

heavy bioturbation 

Silty shale with silt, horizontally 
interbedded to interlaminated, minor 
ripple cross-laminated; heavy 
bioturbation 


Top of McMurray Formation 


Local summary sequence of sedimentary facies 
associations, showing two coarsening-upward 


interlaminated; ripple cross-laminated 
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horizontal, cross-, minor ripple cross- 


horizontal laminae, ripple and bidirec- 
tional cross-laminated; few rip-up clasts 
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The term “shelf” cannot be strictly applied to the siliciclastic sediments in a shallow 
epeiric sea because it implies the existence of a distinct shelf edge. However, little work 
has been done to develop a terminology for shallow siliciclastic environments in epeiric 
seas and the development of this terminology is beyond the scope of this study. The 
terminology of shelf sedimentation is used informally to describe relative water depth and 
is not intended to equate epeiric with pericontinental seas. 

Depositional features of the shallow marine environment are highly variable and 
pooriy understood. There is no one general theory of shelf sedimentation since many 
processes and resulting structures are involved (Anderton, 1976; Harms et al, 1975). 
There are no modern analogues for the sedimentclogy of the Lower Cretaceous Period in 
western North America because there are no large modern epeiric seas. Some 
comparisons can be made between the deposits of modern continental and ancient epeiric 
seas because similar sand bodies are created (Walker, 1979b). However, caution must be 
used in the application of modern to ancient sedimentation since hydraulic conditions and 
basin geometry in continental and epeiric seas are not identical. Furthermore, much of the 
sediment in modern continental seas is reiict glaciofluvial sediment rather than first-cycle 
sediment. Consequently, the relict paleotopographic structures inherited from previous 
sedimentation may affect subsequent sedimentation in modern sediments (Swift, 1969; 
Selley, 1978). 

Several modern and ancient partial analogues (Table 2) have been used in the 
interpretation of the Clearwater Formation in the present study. Partial analogues are 
indirect comparisons of processes that produce similar deposits. 

The shallow marine environment of the Clearwater Sea generated a compiex of 
sedimentary types and grain size distributions, rather than a simple seaward-fining 
sequence. The shelf zone can be considered as transitional between coastal sands and 
deep shelf muds, and the sediments are usually clayey silt to silty sand. 

A number of partially independant factors influence the sedimentary facies 
associations in shelf deposits. These include basin geometry and topography (Mocers, 
1976), rate and type of sediment supply, type and intensity of hydraulic regime, sea level 


flucuations, climate, chemical factors and animal-sediment interactions (Johnson, 1978). 
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Table 2. 


Partial Analogues used in Facies Association Interpretation 


Modern Partial Analogues: 


Location Source 

North Sea Stride, 1963. 
Florida and Bahamas Balt. 1967, 

North Sea Reineck, 1967, 

North Sea Caston, 1972. 
Various - review Ginsberg, 1975. 
Atlantic and North Sea Brenner, 1980. 
Bering Sea Field et al, 1981. 
Washington Continental Sea Nittrouer and Sternberg, 1981. 
New York Bight Vincent et al, 1981. 
North Sea Harms et al, 1982. 


Ancient Partial Analogues: 


Location Geologic Period Source 

Carolina, U.S.A. Upper Cretaceous Swift, 1969. 

Wyoming and Montana Upper Jurassic Brenner and Davies, 1973. 
Various - review Various - review Ginsberg, 1975. 

Wyoming Upper Cretaceous Harms et al, 19/5. 
Scotland PreCambrian Anderton, 1976. 

Wyoming Upper Cretaceous Brenner, 1978. 

Wyoming and Montana Jurassic - Cretaceous Brenner, 1980. 

Norway PreCambrian Levell, 1980. 

Western Interior, N.A. Lower Cretaceous Bridges, 1982. 

Western Interior, U.S.A. Upper Cretaceous Harms et al, 1982. 

Texas Lower Cretaceous Hobday and Morton, 1984. 
Montana Upper Cretaceous Rice, 1984. 

Montana Upper Cretaceous Shurr, 1984, 

Western Interior, U.S.A. Upper Cretaceous Swift and Rice, 1984. 
Wyoming Upper Cretaceous Tillman and Martinsen, 


1984. 
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The hydraulic regime in the shallow marine setting is very complex, consisting of a 
number of interactive currents. Several depositional processes, such as tidal currents, 
meterologic currents (both storm and rip currents}, storm density currents, regional ocean 
currents, and various combinations of the above transport and rework the sediments 
(Howard, 1972; Johnson, 1978; Walker, 197Qb; Harms et al, 1982). 

There is discussion as to the relative significance of storm activity and tidal activity 
in the deposition of sediment in epeiric seas. Storms cause much of the mass movement 
of sediment in shallow marine environments (Brenner and Davies, 1973). Tides strongly 
influence the reworking of these sediments (Off, 1963; Johnson and Belderson, 1969; 
Swift, 1969; Klein, 1977). The tidal range during the Cretaceous Period has been 
estimated at 0.86 m with current velocities up to 10 cm/sec (Slater, 1984). When the 
bottom gradients are very low, tides can affect a large area. The sediments in the 
Clearwater Formation exhibit both storm and tidal depositional features as will be 
discussed in the following section. 

Tectonic rather than eustatic sez level changes also affected sedimentation. The 
epeiric sea was located in a foreland basin adjacent to the developing Western Cordillera. 
This basin subsided in response to the loading of successive thrust plates in the Western 
Cordillera and in response to the weight of sediment being shed from these plates. Each 
thrust created a pulse of subsidence and sedimentation in the foreland basin (Jordon, 
1981; Beaumont, 1981). A net sea level rise of 70 metres during Mannville time is 
predicted due to tectonic activity (Beaumont, 1981). Both orogenic activity and the 
pre-Cretaceous topography control sedimentation in the Clearwater Formation in Cold 
Lake. 

The exact nature of the input mechanism for the Clearwater deposits is not known 
because no paleoshore features were preserved in the area of study. The distribution of 
the sands and the presence of higher amounts of silt and clay-sized material to the 
northwest suggest that the source of sediment was from the southeast. 

Input of sediment may have been from a deltaic complex or from storm surge-ebb 
erosion of a generally prograding shoreline, similar to that described by Banks (1973) in 
the late Pre-Cambrian Innerelv Member of Norway. There is no concrete evidence for a 


well-developed deltaic complex in the Clearwater Formation. There is no evidence of 
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deposition of relatively thick units representing prodelta, delta front and delta top 
sequences. The alternation of sand, silt and shale found in most facies associations 
suggests that fluctuating currents played an important role in deposition and modification 
of the seaiments. The presence of silty shale intraciasts in many of the facies associations 
also suggests periods in which high energy conditions existed. Consequently, it is 
proposed that the sediment supplied to the area was distributed by storm and other 
currents which were strong enough to preciude the development of a delta. The 


coarsening upward sequences are the result of general progradation of the shoreline. 


Coarsening-upward Sequence: Stratigraphic Units J to F 

This coarsening-upward sequence is similar to that described by Harms, et al 
(1982) as a series of broad sandy shoals (stratigraphic unit F) built upward from a floor of 
muddy shelf deposits (stratigraphic unit J). The coarsening-upward trend is seen in the 
geophysical well logs and in the core analysis. This represents an increase in the intensity 
of the currents which were able to carry a progressively coarser bedload. The sediment 
bypassed the nearshore environment iargely due to density stratification currents 
generated by storms. Once deposited on the shelf, the sediment was reworked by 
currents. There was some influence by tidal currents indicatea by widely scattered tidal 
current-generated structures, such as bi-directional low-angle cross-lamination. The 
textural and compositional immaturity of the sands suggests that tidal reworking and 


winnowing was minimal. 


Stratigraphic Unit J: Facies Association 2 

Stratigraphic unit J consists of a bioturbated silty shale with rare moderately 
oil stained argillaceous silt and very fine-grained sand interlaminae. 

Stratigraphic unit J was deposited on the muddy outer shelf. The muds were 
deposited from the suspended load of an river which bypassed the coastal region or 
by storms. The predominately horizontally laminated nature of the sediment is 
suggestive of deposition below effective wave base. Minor tidal and oceanic 
current influences are present in the southeast of the study area and are 
represented by bi-directional low angle cross-lamination and ripple cross-lamination. 


Sedimentary structures, the presence of greater amounts of carbonaceous debris, 
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the coarser-grained silty nature of well 6-3-66-—5W4 and the overall thickening of 
the unit to the northwest indicate proximity to a shoreline towards the southeast. 
The occurrence of traces of Planolites nicholson (1873) is not in itself diagnostic of 
shelf sedimentation, but its presence in all cores suggests a low energy environment 


with a good food supply which supports the interpretation. 


Stratigraphic Unit |: Facies Association 2 

Stratigraphic unit | consists of a bioturbated silty shale with interbeds of silt 
and very fine-grained sand. 

Sedimentation in this shelf facies was dominated by storm-generated density 
stratified currents. The intermittent storms divided the sediment into bedload and 
suspended load and moved it across the shelf. Similar sequences have been 
interpreted by Walker (1979b) for the Cardium Formation of Alberta and Swift and 
Rice (1984) for the Upper Cretaceous of Montana and Wyoming. The coarse 
sediments represent a periodic short-lived event when the silt and very fine-grained 
sand was introduced into an otherwise quiescent environment. Successive storms 
built up a silt and mud sequence. Between storms the silty shales were consolidated 
and bioturbated. This accounts for the relative paucity of storm rip-up clasts in the 
storm deposited sands. The high degree of bioturbation suggests that storm events 
were sufficiently far apart to allow for the establishment of large communities of 
organisms which fed on these deposits. After each storm event the communities | 
re-established themselves. The more seaward portion of the stratigraphic unit is 


situated to the northwest where it is more shaly. 


Stratigraphic Unit H: Facies Association 3 

Stratigraphic unit H consists of a very fine-grained sand with interlaminae to 
interbeds of silty shale and silt. The silty shale is heavily bioturbated at the base of 
stratigraphic unit H and decreases upward. 

The increased strength of storm currents in a shallow water environment is 
indicated by the sedimentation pattern of stratigraphic unit H. The density stratified 
storm currents were capable of transporting larger amounts of very fine-grained 


sand as bedload. The greater number of silty shale intraciasts in the sands suggests 
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that the silty shale did not deveiop enough cohesion to withstand the next storm 
pulse as effectively as in stratigraphic unit |. The fewer number of traces of 
Planolites nichoison (1873) suggest that higher energy conditions did not allow the 
establishment of large communities of organisms. This interpretation is reinforced 
by the presence of load casts indicating differential loading of the sands on 
water-saturated muds during deposition. 

The presence of the trace of Skolithos linearis haldeman (1840) is not 
paleoenvironmentally diagnostic. Although many researchers (Seilacher, 1967; 
Crimes, 1970; Frey, 1975) claim that the presence of the trace of this suspension 
feeder indicates a nearshore environment, others (Alpert, 1974; Basan and Frey, 
1977) point out that this trace can be found in environments ranging from the outer 
shelf to the shoreface. These traces were made by organisms particularly adapted 
to unstable environments which are characterized by rapid sedimentation. The 


association of Skolithos linearis haldeman (1840) and Planolites nicholson (1873) is 


suggestive of an oxygen-depleted stressed environment (Ekdale et al, 1984). 


Stratigraphic unit G: Facies Association 3 

Stratigraphic unit G consists of very fine-grained sand with interlaminae to 
interbeds of silty shale and silt. Sand units are thicker than those in stratigraphic unit 
H. Bioturbation is less extensive than in stratigraphic unit H and decreases upward. 

The structures and lithology of this stratigraphic unit suggest that it was 
deposited closer to the paleoshore in a shelf setting. The silty shale and sand 
interbeds are suggestive of storm surge-ebb deposition. Sand was deposited during 
the waning period of a storm and the silty shale was deposited and partially 
bioturbated during quiescent periods. The isopach map pattern (Figure 18) suggests 
that the pre-Cretaceous highs served as loci for stratigraphic unit G deposition. The 
abundance of organic material, especially in well 6-3-66-5W4, is suggestive of 
proximity to the paleoshore to the southeast. Some tidal current influence is 
indicated by the presence of bi-directional low angle cross-laminae. Inclined 
low-angie cross-laminae (which may be trough cross-laminated) suggest current 
reworking of the sediment. An exclusively tidai interpretation for this facies is 


difficult because the tidal-generated structures are found in association with 
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structures which may have been created by oceanic currents or waves. 
Biogenic activity is sparse, indicating a slightly higher energy level and/or a 


more rapidly prograding substrate. 


Stratigraphic Unit F: Facies Associations 4 and 6 

Stratigraphic subunit F, (Facies Association 6) consists of finely laminated 
poorly sorted sand with rare carbonaceous silty shale laminae. Stratigraphic subunit 
F, (Facies Association 4) consists of a very fine-grained, moderately poorly sorted, 
laminated to inclined low-angle cross-laminated sand unit with silty shale and 
siltstone interbeds. Stratigraphic subunit F, (Facies Association 6) consists of 
horizontally laminated very fine-grained, poorly sorted sand with rare silty shale 
interlaminae. 

Stratigraphic unit F represents large amounts of storm-transported very 
fine-grained sand, silt and clay modified by tidal, oceanic and subsequent storm 
currents. Similar sequences have been reported by Ball (1967) in recent sediment in 
Florida and the Bahamas, Amajor (1980) in the Lower Cretaceous Viking Sandstone 
of central Alberta, Brenner (1980) in the Jurassic and Cretaceous of Wyoming and 
Montana, Field et al (1981) in the Bering Sea, Vincent et al (1981) in the New York 
Bight, and Swift and Rice (1984) in the Upper Cretaceous of the western interior of 
the United States. Some of the sediment may have been transported during peak 
tides as well. 

The coarsening-upward nature of the deposit and the thicker accumulation of 
sand may be due in part to shallowing as the inferred shoreline located towards the 
southeast prograded. The current-generated structures, the interlaminations of silty 
shales and the lack of bioturbation suggest that higher energy conditions existed at 
the time of deposition of stratigraphic unit F. Thick (up to 12 m) units of inclined 
low-angle cross-laminated sands suggest the formation and maintenance of large 
scale structures. The isopach map (Figure 19) shows a northwest to southeast trend 
to sediment distribution. This suggests that after initial deposition, the sediment was 
reworked and loca! topographically higher structures were created and maintained. it 
is proposed that stratigraphic unit F was deposited in the form of poorly developed 


bars and ridges. 
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The classification of sand deposits in the shallow marine environment is 
confused (Levell, 1980; Belderson et al, 1982) and still in debate. The classification 
proposed by Shurr (1984) is based primarily on sand body geometry. According to 
this classification, the deposits in the Clearwater Formation are sand bars to sand 
ridges in well 6-5—66—-5W4 where there is an accumulation of 12 m (40 ft) of 
cross-laminated to trough cross-laminated sand. Additional evidence is supplied by 
the spontaneous potential and gamma ray responses (Serra and Sulpice, 1975; 
Pirson, 198 1b) both of which are subdued. This suggests the presence of bar-type 
deposits. 

The deposits of stratigraphic unit F are similar to those described by Rice 
(1984) and Tillman and Martinsen (1984) in the Upper Cretaceous Shannon 
Sandstone of the northwestern plains area in the United States. The deposits were 
described as elongate discontinuous sand bodies on a shallow flat shelf with 
associated bar margin and interbar deposits. The deposits in stratigraphic unit F are 
also similar to those described by Brenner (1978) in the Sussex Sandstone of 
Wyoming as a shelf ridge complex. These deposits aré comparable to the storm and 
tide built sand ridges on the modern Atiantic Continental Shelf or the Bight of the 
North Sea. However, the deposits found in Cretaceous sediments are on a smalier 
scale than the modern ridge because the epeiric sea was much more shallow which 
limited the vertical growth of the ridges. The presence of the coarsening-upward 
ridges reflects a shallowing of the Lower Cretaceous epeiric sea and an increase in 
the supply of silisiclastic material to the area. The occurrence of these bars and 
ridges was controlled by topographic highs on the sea bottom. Their orientation was 
controlled by southeast-northwest trending currents (Figure 19). Slight topographic 
variations can result in spatial velocity gradients that induce deposition. Once 
deposition was initiated, the bars enhanced the topography and became 
self-perpetuating closed systems. 

Thinner accumulations of very fine-grained sand and silty shale represent 
smaller scale sand waves and ripples superimposed on the bars and ridges. The 
structures exhibited include ripple cross-lamination, clay drapes on ripples and 


bi-directional low-angle cross-lamination. These sedimentary structures reflect 
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some tidal reworking of the sediment, although it is difficult to distinguish tidal 
current from other current deposits. Flaser structures found in this facies can be 
interpreted as either storm-generated or tidal-generated. 

The interbedded finely laminated and ripple cross-laminated silty shale, 
argillaceous silt and very fine-grained sand found in stratigraphic unit F represent 
interbar and bar margin deposition. 

Erosional surfaces represent uncorrelative channels. The presence of such 
channels is suggestive of proximity to the paleoshore. These channels are infilled by 
ripple cross-laminated and dune cross-laminated sands which migrated on the 
channel floor. 

There is no evidence of emergence in these deposits, although wells 
6-3-66-5W4 and 6-5-66-5W4 contain larger amounts of carbonaceous debris, 


indicating their proximity to the paleoshore. 


Coarsening-upward Sequence: Stratigraphic Units E to B 

This coarsening-upward sequence represents deposition in an inner shelf 
environment. The inner shelf environment is characterized by a progression towards a 
coarser sand-size fraction, a decrease in bioturbation and an increase in sorting of the 
sediment due to an increase in energy shoreward. 

The inner shelf is dominated by wind (storm) effects unless there is a high tidal 
influence (Mooers, 1976; Harms etal, 1982). 

No foreshore or upper shoreface sedimentation was recognized in the cores and a 
paleoshore trend cannot be reconstructed. There are several explanations for this 
omission in the rock record. In this shallow marine environment with its low-sloping sea 
bottom any flucuation in sea level will rapidly shift the shoreline over a great distance. This 
would preclude the development of distinct beach sequences. The effects of storms and 
tides would cause erosion of any beach deposits. This creates a paucity of shoreline 
deposits and a trend towards sand storage in submarine sand bodies rather than beach 
deposits. The zone of waves and tides in an epeiric sea may have been miles from the 


shore, preventing the the formation of beach deposits. 
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Stratigraphic Unit E: Facies Association 6 

The basal portion of stratigraphic unit E consists of horizontal, finely 
interlaminated sand and silt. The upper portion consists of horizontally laminated and 
ripple cross-laminated silty shale and very fine-grained sand. 

The deposit in stratigraphic unit E is the best indicator of tidal influence in the 
Clearwater Formation. The rhythmic finely interlaminated sediments are tidal bundles 
reflecting small scale repeated alteration in sediment transport conditions. The sand 
laminae were deposited as bedload during the higher velocity phase of the tidal cycle 
and are the remnants of very flat-crested rippies. The silt laminae were deposited 
during the slack phase of the tidal cycle when fines came out of suspension. 

The presence of tidal bundles is not paleoenvironmentally diagnostic within 
tidal settings as they can be found in both intertidal and subtidal environments. The 
lack of emergent structures, such as mud cracks, suggests that the tidal bundles 
were deposited in a subtidal environment. 

The presence of less well-developed tidal bundles and bi-directional 
low-angle cross-laminae in well 2—-19-66-5W4 suggest that the effects of tidal 
currents were widely distributed. However, the lack of deposition of tidal bundles in 
wells 10-7-66-5W4 and 10-—17-66-—5W4 and the inability to trace this 
stratigraphic unit regionally suggest that tidal currents of sufficient regularity to 


produce tidal bundles were aerially restricted phenomena. 


Stratigraphic Unit D: Facies Associations 1 and 3 

Stratigraphic subunit D, (Facies Association 1) consists of a heavily 
bioturbated silty shale. Stratigraphic subunit D, (Facies Association 3) consists of a 
bioturbated poorly sorted sand with interbeds of bioturbated silty shale. 
Stratigraphic subunit D, (Facies Association 1) consists of a heavily bioturbated silty 
shale. 

A change in sediment supply is indicated by the deposition of stratigraphic 
unit D. The thick accumulations of bioturbated silty shale suggest that the supply of 
very fine-grained sand was interrupted. Mud zones can have a number of 
interpretations, such as a muddy coast, a nearshore mud belt, an mid-shelf mud belt, 


an outer-sheif mud belt or a mud blanket off a delta (McCave, 1972; Anderton, 
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1S76; Johnson, 1978). With the exception of the outer-shelf mud belt, all of these 
features are of limited aerial extent. In the case of stratigraphic unit D, the shales 
have a broad distribution beyond the area of the regional study (Wightman, personal 
communication). The isopach map (Figure 22) suggests that thicker accumulations 
occur seaward to the north and west. This would suggest that broader forces, such 
as regional subsidence and/or a temporary decrease in coarser sediment supply due 
to a temporary cessation of tectonic activity in the Western Cordillera, which was 
the main source of detritus, were in effect at this time. This caused the development 
of an extensive shallow marine shelf mud. 

The presence of intraclasts suggests that the interbed of very fine-grained 
sand represents storm-transported sediment. However, the highiy bioturbated 
nature of the sana and the presence of the silty shale above it suggest that after 
sand deposition, relatively quiescent conditions continued. This allowed the 


extensive biogenic reworking of these sediments. 


Stratigraphic Unit C: Facies Association 5 

Stratigraphic unit C consists of graded sand, silt and shale interbeds. 

The deposit in stratigraphic unit C represents a return to the previous very 
fine-grained sand, silt and clay sediments derived from a source situated to the 
southeast. 

Stratigraphic unit C is interpreted as being located in an inner shelf 
environment between the shoreface and deep water. Proximity to the paleoshore is 
suggested by the high organic content in the deposit, although there is no evidence 
of emergence. Storm and tidal influences were dominant in this environment. 
Storm-generated turbidity current deposits were common and were separated from 
each other by finer-grained deposits and wave ripple cross-laminated sands. 
Periodic storm turbidity currents transported the very fine-grained sand, silty shale 
intraclasts, silts and clays as suspended load creating a graded sand, silt and shale 
sequence of up to 0.37 m (1.25 ft) in thickness displaying Bouma sequences. The 
sediment load of the Western Interior Cretaceous seaway was very fine-grained anc 


created sequences skewed towards the fine-grained end. 
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The deposits in stratigraphic unit C also reflect deposition from tidal 
currents, exhibited by bi-directional low-angle cross-lamination which is found in all 
cores studied. The inclined low-angle interlaminated silty shale, silt and sand 
represent accretion surfaces of a small scale bedform. 

Weak bioturbation is indicative of fast migration of sediments in an inner 


shelf environment. 


Stratigraphic Unit B: Facies Association 6 

Stratigraphic unit B consists of poorly sorted laminated sand interlaminated 
to interbedded with laminated silty shale. 

The accumulation of up to 2.1 m (7 ft) of horizontally laminated to inclined 
cross-laminated sands in stratigraphic unit B represents the development of small 
scale shoaling sand bars above the effective wave base in an inner shelf 
environment. The isopach map (Figure 24) outlines these bars and suggests that they 
were oriented parallel to the inferred shoreline to the southeast. This feature is best 
developed in well 10-—17—-66—-5W4. The interlaminated silty shale represents bar 
margin and interbar sedimentation in the lee of these bars. Shallow marine tidal or 
storm-generated dunes are represented by the thinner accumulations of 
cross-laminated sands associated with laminated sand and silty shale. 

There is no evidence of emergence in these deposits, although large 
concentrations of carbonaceous debris indicate proximity to the paleoshore. Weak 
bioturbation is indicative of the shallow water conditions and fast migration of 


sediments. 


Stratigraphic Unit A 

This stratigraphic unit was not cored so definite paleoenvironmental 
interpretation cannot be made. Equivalent shales in the Cold Lake, Wabasca and 
Athabasca oil sands deposits have been interpreted as open marine outer shelf muds 
(Miossop et al, 1981 and others). This stratigraphic unit represents a transgression 


of the epeiric sea over the inner shelf sediments. 
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Summary 

The Clearwater Formation in the study area consists of two coarsening-upward 
sequences (Figure 26). The first sequence (stratigraphic unit J to stratigraphic unit F) was 
deposited in an outer shelf environment and the second sequence (stratigraphic unit E to 
stratigraphic unit B) was deposited in an inner shelf environment. 

The first coarsening-upward sequence is approximately three times thicker than 
the second coarsening-upward sequence. This apparent discrepancy in thickness can be 
attributed to the presence of stacked offshore bars in stratigraphic unit F. The silty shale 
and siltstone in stratigraphic subunit F, represent interbar or bar margin deposits which 
delineate the sand bars. The inner shelf sediments (stratigraphic units B and C) were 
deposited after a period of subsidence. This subsidence resulted in a change in the basin 
geometry and the sediment source area Consequently, the two coarsening-upward 
sequences represent different sets of processes and products and do not reflect a single 


prograding shoreline. 
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IV. Mineralogy 


A. Bulk Mineralogy 

Forty-one thin sections in well 10-17-66-5W4 were examined to establish the 
mineralogy of the Clearwater Formation at this location. Relative percentages of point 
counts are given in Appendix ||. Mineralogic composition is summarized and compared 
with Putnam and Pedskalny (1982) and Harrison et al (1981) in Table 3. 

Grain mineralogy consists of quartz with lesser amounts of chert and feldspar. 
Igneous, volcanic, metamorphic and sedimentary rock fragments were present in most 
samples. The complex mineralogy suggests a combination of source rocks. The rock is 
composed of feldspathic litharenites to litharenites, according to Folk’s (I968) 
classification (Figure 27). Points plotting low on the feldspar / igneous rock fragments end 
of the ternary diagram contain a calcarous matrix which has replaced many feldspar grains 
as well as other primary mineralogy. Ee 

Composition compares well with Harrison et al’s (1981) study in Townships 64 and 
65, Range 3W4 and Putnam and Pedskalny’s (1983) study in Township 64, Range 4W4, 
both of which are situated southeast of the present study. However, the B.P. pilot plant 
site exhibits a more equitable distribution of feldspathic litharenites and litharenites and 
contains slighty more quartz grains. This is due to farther transport of the sediments to the 
B.P. pilot plant site. While there are no obvious trends in the distribution of the two 
compositions, the upper stratigraphic units (B, C and D,) and the lower stratigraphic units 
(G, H, and |) are composed predominantly of feldspathic litharenites. In stratigraphic units E 
and F, litharenites predominate, although feldspathic litharenites are a component of 
stratigraphic subunits F, and F,. 

The dominant lithic grains are metamorphic rock fragments (29.5 percent). This 
finding is at variance with Harrison et al (1981) and Putnam and Pedskalny (1983) who 
found that volcanic rock fragments and chert were the dominant lithic grains. The main 
reason for this discrepancy is that the large proportion of polycrystalline quartz rock 
fragments were classified as volcanic rock fragments by Putnam and Pedskalny (1983) and 
presumably by Harrison et a! (1981), while the present study classifies them as 


metamorphic rock fragments. This discrepancy will be discussed in detail under 
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Tables; COMPARATIVE MINERALOGY 


THIS STUDY 
Putnam 
Overall Recalculated and Harrison 
Composition to 100% Pedskalny et al 
(%) (1982) (1981) 
(%) (2) 

Quartz here VA S323 2o08 2120 
Chert 4.0 Se) bie4 20.0 
Feldspar 4.7 FO.0 27.9 20.0 
Igneous Rock Fragments 3.4 Tom = = 
Metamorphic Rock Fragments 1339 2975 10.0 se 5) 
Sedimentary Rock Fragments raat 4.5 11.4 ced 
Volcanic Rock Fragments GE: 14700 Jie) asi 
Clastic Carbonates 4.0 #0030 UES 
Glauconite 253 
Mica 05 
Chlorite Zao 
Sericite G25 
Zircon 0.8 
Opaques 0.5 
Organics 7.4 
Clay 4.0 
0i1 1E2 
Pore Space 4.4 
Matrix ae & 
Calcareous Cement 9./ 
Quartz Overgrowths Ue 
Unidentified Rock Fragments O22 
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QUARTZ 


ee METAQUARTZITE 


SUBLITHARENITE 


ARK OSE LITHIC ARKOSE FELDSPATHIC 
LITHARENITE 


LITHARENITE 


o Calcareous cemented samples 


e Non-cemented samples 


FIGURE 27. Ternary plot of point-counted thin sections from 
well 10-17-66-5W4, showing classification of sands 
as feldspathic litharenites and litharenites. 


(after Folk, 1968) 
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metamorphic rock fragments. Note that the classification of polycrystalline quartz rock 
fragments does not affect the ternary diagram, since both metamorphic and volcanic rock 


fragments plot as labile rock fragments. 


Quartz 

Unicrystalline quartz grains make up between | percent and 24 percent of the tota! 
mineralogy. They are generally subangular to subrounded and uniformly very fine-grained. 
Some very well-rounded grains were observed in the massive sands of stratigraphic unit 
F, suggesting reworking of this sediment. Straight and undulous or strained extinction was 
observed in the quartz grains. 

Authigenic quartz overgrowths (Plate 6-2) are a minor constituent, but were found 
in most thin sections examined. Most of these overgrowths exhibit weathering and are 
derived from a previous depositional cycle. 

Authigenic feldspar overgrowths (Plate 4-1) are rare and show abrasion, indicating 


formation in a previous depositional cycle. 


Chert 

Chert grains make up between | percent and 7 percent of total rock volume. This 
can only be considered an estimate, since it is difficult to differentiate between some 
polycrystalline quartz grains and chert. An arbitrary distinction was made which assigned 
grains with randomly oriented microcrystals less than 10 microns in diameter to chert 
(Folk, 1968, p. 70). The colour of chert ranges from white to black. Most grains are 
subangular, although rare rounded grains and lathe-shaped fragments were noted in 


stratigraphic unit B. 


Feldspar 

Feldspars, both plagioclase and potassium, are found throughout the core. They 
make up between | percent and 9 percent of tota! rock volume. Many of the plagioclase 
feldspars exhibit a euhedral crystal outline, well-defined crystal zoning and albite twinning 
(Plate 4—2), all indicative of a volcanic origin (Scholle, 1979). Potassium feldspars are 


dominantly orthoclase. Most feldspar grains are subangular with some platy fragments 
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noted. This suggests moderate transport distances. Overburden pressure is sufficient to 
fracture some feldspar grains (Plate 4—3). Alteration of the feldspars to kaolinte clay and 
calcite is seen, but many feldspar fragments are quite fresh, a reflection of little 


reworking and gocd preservation of the sediment. 


Rock Fragments 

Plutonic (igneous and volcanic), metamorphic and sedimentary rock fragments are 
found throughout the core and can contribute up to 38 percent of the total rock volume. 

Igneous rock fragments were distinguished from metamorphic rock fragments on 
the bases of the presence of igneous trachytic texture, coarse equant crystals and unequal 
crystal sizes in the grain (Scholle, 1979). Many grains are highly corroded. The proportion 
of igneous rock fragments increases from none at the base of the core to 13 percent in 
stratigraphic unit B which is situated near the top of the formation. This suggests a shift in 
source area during Clearwater Formation time. 

While it is agreed that the immature nature of the rock composition indicates that 
the major source of detritus is from the newly-forming Cordillera to the west (Williams, 
1963; Maycock, 1964; Rudkin, 1964; Cameron, 1965; Rapson, 1965; Putnam and 
Pedskalny, 1983), the presence of the degraded igneous fragments may suggest a minor 
component of detritus from the Shield area to the northeast (Mellon, 1967; Keeler, 1980; 
Harrison et al, 1981). Further petrographic work is needed to confirm this proposition. 

Volcanic rock fragments, derived from the Cordilleran region to the west (Putnam 
and Pedskalny, 1983) are found throughout the core in relative percentages ranging from 
1 percent to 8 percent. Searlesite fragments were observed in stratigraphic unit B. The 
identification of volcanic rock fragments was based primarily on the presence of 
predominately feldspathic micophenocrysts. These microphenocrysts are commonly 
untwinned and strongly zoned (Plate 4—4)}. 

Metamorphic rock fragments contribute between 3 percent and 26 percent to the 
mineralogy. Grains exhibit strongly oriented, sutured and rare schistose textures and are 
subangular to subrounded, indicative of moderate transport. Polycrystalline quartz was 
identified as a metamorphic rock fragment (Plate 5-1) rather than a volcanic rock fragment 


(c.f., Putnam and Pedskalny, 1983). Quartz from a volcanic source is usually 
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monocrystalline and fractured (Blatt et al, 1980, p. 289). Rapson (1965) could not 
recognize volcanic quartz in age equivalent sediments deposited in the southern Rocky 
Mountains of Alberta. It is proposed, therefore, that polycrystalline or composite quartz is 
indicative of a metamorphic source (Mellon, 1967; Scholle, 1979; Blatt, 1982). 
Metamorphic quartz exhibiting inclusions was also noted (Plate 7-3). 

Fine-grained sedimentary rock fragments of silt, silty shale, and shale contribute 
between 1 percent and 5 percent to total rock volume. Most grains exhibit microscopic 
orientation, reflecting original lamination. Another distinguishing characteristic is plastic 
deformation due to compaction (Plate 7-1). 

Clastic carbonates are found throughout the core. These grains contribute 
between 1 percent and 5 percent to total rock volume and are subangular to angular, 
indicating short transport distances. In a4 m (12 ft) long portion of core in stratigraphic 
unit E and the upper portion of stratigraphic unit F, clastic carbonate concentrations of up 
to 24 percent are found. Increased concentrations have been reported from other wells in 
the area (Wightman, personal communication). The clastic carbonates are derived from 
erosion of pre-Cretaceous highs in the area. The high concentrations are probably 
reflective of increased tectonic activity resulting in uplift of the pre-Cretaceous islands 
which then shed more carbonate detritus into the basin at the time of deposition of 


stratigraphic unit E and the upper portion of stratigraphic unit F. 


Glauconite 

Glauconite is common throughout the core and can attain concentrations of up to 
10 percent of total rock volume, although concentrations average 4 percent. Most pellets 
are well rounded, slightly larger than surrounding grains and occasionally deformed due to 
compaction. Pellet morphology is that of singular spheroid or ovoid grains (Plate 5-3). The 
pellets have the structure of an illite-smectite clay and can be derived from these clays or 
micas, mud fillings of organisms, or fecal pellets (Burst, 1956; Triplehorn, 1966). 
Glauconite occurs in both detrital (Johnson, 1978) and early authigenic forms, replacing 
sedimentary rock fragments and biogenic pellets. It is formed during and shortiy after 
deposition (Mellon, 1967) at the sediment-water interface (Wilson and Pitman, 1977). An 


example of glauconite replacing original oolites can be seen in Piate 5-3. Most glauconite 
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found in the Clearwater Formation appears to have been formed in situ. Its presence 


reinforces a marine origin for these sediments. 


Calcareous Cement 

Calcareous cement is a diagenetic product found in discrete beds and lenses. It has 
several habits. The cement is found replacing original matrix and most mineral grains (Plate 
5-4) or as calcareous veins displacing the original rock framework (Plate 6—1). Although 
this cannot be seen in core, it is assumed that the calcareous cement is nodular. At best, it 
can be correlated over short distances as elongate lenses. 

Cone-in-cone structure is found in stratigraphic subunit F,, in both wells 
2-19-66-5W4 and 10-17-66-5W4. This structure is described as an early diagenetic 
product resulting from differential pressure during crystal growth (Fairbridge and 
Bourgeois, 1978). It is characterized by an abundance of tight circular cones which stand 
with the cone axis perpendicular to the cone-in-cone layer (Pettijohn, 1975). 

The structures and occurrence of cone-in-cone in the Clearwater Formation are 
similar to those described by Gilman and Metzer (19687) from the Upper Devonian 
Canadaway Formation. In both cases there are two sets of cones, one pointing upward and 
the other downward, separated by a layer of calcareous cemented siltstone. Fractures 
infilled with calcareous cement are associated with the cone-in-cone structures. Gilman 
and Metzer (1967) postulate that these structures originate from a syngenetic concretion 
of fibrous aragonite caused by pressure during dewatering of the surrounding sediments. 

Siderite cemented nodules with diameters of up to 4 cm occur rarely in 
stratigraphic units | and J. This is indicative of reducing conditions during formation of 


these nodules (Scholle, 1979). 


Organic Debris 

Organic fragments are ubiquitous, but exhibit a highly variable distribution, ranging 
up to 33 percent at the 442.9 m (1452 ft) depth of stratigraphic unit C and as low as 2 
percent in the basal portion of the same stratigraphic unit. Several! small coal laminae were 
observed in core in basal stratigraphic unit B and upper stratigraphic unit C. A similar 


variable distribution can be seen in most other stratigraphic units. Organic fragments and 
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laminae show a preferred orientation in thin section (Plate 6-1) and are commonly 
distorted due to compaction. The organic fragments are usually much larger than 
surrounding mineral grains (Plate 6—2) due to their lower specific gravity. 

Presence of organics usually reflects proximity to a shoreline, but due to their low 
specific gravity, organics can be transported far out onto the shelf as well. In stratigraphic 
units B and C, the high organic content may be indicative of promixity to a shoreline or 
relatively higher organic input at this time. In other stratigraphic units, the organic 


component is explained as normal shelf input in an area of high organic activity on land. 


Accessories 

Small amounts of both biotite and muscovite mica (Plate 7— 1), hematite and chlorite 
grains (Plate 5—2) were observed in thin section. Relative percentages never exceed 2 
percent. Chlorite grains exhibit both ultra-blue birefringence and normal low birefringence 
(Scholle, 1979). Very fine-grained zircon is the most abundant non-opaque heavy mineral 
in this area, ranging up to 6 percent in stratigraphic unit F. Dead oil is also quite common as 
a pore lining and as a stain on clay (Plate 7—3). Differentiation of oil and clay is difficult in 
thin section (Scholle, 1979). The authigenic clay was formed subsequent to burial of the 
sediments by precipitation from fluids or the alteration of unstable detrital grains. 
Examples of feldspar partially altered to kaolinite and biotite mica partially altered to 
chlorite were seen in thin section. 

Authigenic pyrite is commonly seen as discrete cubes lining pore spaces along with 
clays. It is also found as irregular and frambroidal masses, replacing grains and organic 


matter (Plate 7-2). 


B. Clay Mineralogy 

The Clearwater Formation in well 10-—17-—66—5W4 is composed of a complex 
suite of detrital and authigenic clays. 

Clay mineralogy was determined by X-ray diffraction for thirty-one oil sands and 
five shale samples (Table 4). Crushed samples were centrifuged to obtain a less than 2 
micrometer fraction and smear slides were made. Humidity conditions were controlled by 


treating the samples with ethylene glycol. The samples were then analyzed using a 
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Table 4. Relative Percentages of Clay Minerals by X-ray Diffraction 
(<2um size fraction) 
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Chlorite 
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diffractometer. When the presence of chlorite in the clay mineral assemblage was 
suspected, the sample was heated and analyzed by the diffractometer again. The heating 
treatment separates the kaolinite and chlorite peaks. 

Relative percentages of kaolinite, illite, smectite and chlorite were determined by 
the peak intensity method. This method gives relative percentages comparable with those 
generated by the weighted peak-area method also in use (Scafe, personal communication). 

Detailed Scanning Electron Microscopy to determine authigenesis of the clays was 
beyond the scope of this study. 

Relative percentage clay types are highly variable, but the distribution of clay types 
relates well to the lithology of the samples. Sand, argillaceous sand, mixed sand and clay, 


and shale samples can be differentiated on the basis of their clay content. 


Sand Samples 

Oil sands samples are characterized by a very fine, grain-supported texture with 
good oil saturation. Samples never exceed four percent clay by sieve analysis. The clays 
are primarily authigenic (Wilson and Pitmann, 1977). With the exception of samples 
80-S-11 and 80-39-16, all samples contain a high proportion of kaolinite. Illite is found in 
low amounts except in stratigraphic units D and E and stratigraphic subunit F,, where it is 
found in amounts equivalent to those of kaolinite. Smectite and chlorite are found in lesser 
amounts, but are slightly more significant in the sands of stratigraphic units D and E, 
stratigraphic subunit F, and uppermost stratigraphic subunit F;. 

Some clay was observed in thin section as discrete particles, oil stained booklets 
of kaolinite lining pore spaces and bridges (usually illite) between grains (Neasham, 1977). 
Care was taken to distinguish between genuine bridging and apparent bridging due to 


plucking of grains during thin sectioning. 


Argillaceous Sand Samples 

These samples, containing 1 to 6 percent clay by sieve analysis, occur in close 
proximity with shales or shaie intraclasts. The clays exhibit a tendency towards equal 
amonunts of kaolinite, illite and smectite, although kaolinite still predominates. Chlorite is 


an accessory Clay mineral here. 
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Mixed Sand and Clay Samples 

Matrix clay samples contain comparatively large amounts of silt-sized detrital clays 
(six to thirty-nine percent by sieve analysis). These clays are derived from the crushing and 
disintegration of argillaceous rock fragments (Rapson, 1965). Samples which contained 
large lenses of silt-sized material or were entirely composed of matrix-supported grains 
were placed in this category. 

Due to their silt size and compaction of the sediment, these clays were not 
identifiable in thin section. X-ray diffraction results show that illite and smectite are the 
major constituents, while kaolinite and chlorite are accessory clay minerals. There is an 
association between illite and chlorite in that both seem to increase and decrease at the 


same rate. 


Shale Samples 
The clay in shale interbeds is predominately illite and smectite with chlorite as a 


minor constituent. Kaolinite is an accessory clay mineral in these samples. 
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V. Geophysical Well Log Response in Shaly Sands 


A. Previous Work 

The composition of the bitumen in the Cold Lake oil sands area is considered 
transitional between that of the Athabasca and Lloydminster areas. It has an API gravity of 
between 10° and 14° (Mellon, 1967). Reserve estimates for the 3,500 square mile reserve 
were originally thought to be 75 billion barrels (Mellon, 1967), but have been raised to 164 
billion barrels (Energy Resources Conservation Board, 1973). 

The Clearwater Formation contains the most continuous reservoirs in Cold Lake 
(Vigrass, 1968; Buckles, 1979). The reservoirs here are between 10 mand 15 m thick 
(Matheny, 1979; Mossop et al, 1981) with local thicknesses of up to 35 m (Mossop et al, 
1981). Porosity of the reservoir sands is between 28 and 43 percent, averaging 37 
percent (Jardine, 1974; Matheny, 1979), while permeability is between 300 and 1300 
millidarcies (Jardine, 1974). Oil saturation is 14 to 16 weight percent (Mossop et al, 

1981). Control on the oil saturation has been considered to be primary porosity and 
permeability since post-depositional changes were considered minimal (Mossop et al, 
1981). Several mechanisms and combinations of mechanisms have been proposed to 
account for the accumulations of oil in the Clearwater Formation at Cold Lake. 

Structure, specifically the influence of Paleozoic highs, is one of the major factors 
controlling oil accumulation (Wickenden, 1948; Mellon, 1967; Energy Resources 
Conservation Board, 1973; Carrigy and Kramers, 1974; Jardine, 1974; Mossop et al, 
1981). Salt collapse during and after Mannville time (Jardine, 1974) has also been envoked 
as a structural control on bitumen accumulation. 

In the Clearwater Formation at Cold Lake, major oil occurrences have been located 
on large, structurally high features which coincide with Paleozoic highs, while smal! 30 
foot oil accumulations were found scattered on structural highs downdip to the southwest 
of the major accumulations (Vigrass, 1968). Water is located on the flanks of these 
features (Kendall, 1977). On the eastern edge of the oil sands deposit, synclinal features 
create a gentle eastward dip which marks the termination of oi! accumulation (Energy 


Resources Conservation Board, 1973). 
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Stratigraphy has been identified as another major controi on oil accumulation in the 
Cold Lake oil sands area (Wickenden, 1948; Vigrass, 1968; Carrigy and Kramers, 1974: 
Jardine, 1974; Mossop et al, 1981). The Clearwater Formation, which contains 40 
percent of Cold Lake hydrocarbons (Buckles, 1979), has the thickest, most laterally 
continuous sands (Kendall, 1977; Buckles, 1979; Mossop et al, 1981). However, the 
vertical continuity of the Clearwater Formation is disrupted by shale interbeds and 
cemented siltstones (Kendall, 1977). The clay and shale interbeds and differential 
compaction of sediments (Jardine, 1974) reduce the effectiveness of these reservoirs. 

Petrography also influences oil accumulation (Buckles, 1979; Harrison et al, 1981). 
The complex mineralogy of the deposit influences quantitative log interpretation (Buckles, 
1979). Mean grain size, the amount of fines, and the presence of diagenetic clay all affect 
the emplacement and recovery of the oil (Harrison et al, 1981). The diagenetic clays have 


detrimental effects on the recovery of oil from the oil sands (Gallup, 1974). 


B. Introduction to Geophysical Well Log Study 

The geophysical well log study of the Clearwater Formation highlights the 
problems associated with the interpretation of reservoir characteristics in a shaly sand 
deposit. 

A shaly sand is defined as a sand which contains greater than five percent clay 
(Fertl, 1972). The ultimate aim is to deal with the various clay types, the distribution of the 
clay in the sediment and the effect of the clays on the geophysical parameters measured 
(Helander, 1980). 

Core-derived parameters, such as porosity, weight percent bitumen and weight 
percent water are inaccurate for several reasons which will be discussed in detail below. It 
is difficult if not impossible to obtain accurate values for these parameters from 
unconsolidated sands under conditions truly representive of in situ conditions (Collins, 
1977). Parameters derived from individual well logs are also inaccurate due to the various 
influences of the clays and other minerals on the log response. Consequently, multiple 
porosity log analysis techniques have been developed to reconcile log-derived parameters 
with one another and with those derived through routine core analysis (Tixier et al, 1968; 


Poupon and Gaymard, 1970; Hesiop, 1975). 


Ri ri ris lures he pr ai wi , i ; . ar 
7 
PRR z sell oe a aah 


OD 2ruewsrind Tt oe: i hd hula tei M ‘ f A 
yilgxgial Taoey teeatoxtt ite wai) ies 1 afeios 
ar? “eveyV ci (7 8R7. nace" roTé a ow pee 
’ ‘vie eit . , 
we enachelr sare “ef pei wee ed 1a wands 7 hon es 


_ 
= 7 OL 
init ane ib bre sear! em oi a ‘ gat acer abe 28 Le bere 
vy * ; re an Gi at _ Ww 
sort 7o a2onaviigette 4% oogimt gas anne! _—e ae To fo inex ne 
’ i . i - . a 


vert evel eebtai id) Cer RAS becghiatieuttnl ost aa ae 


78GF js to MORN 
per EW ricita re cy meden, gat oxo Hash? scree ne arto vooiere 
rove Hila vain sitanagait te \ cana ¢ cate Eire : an pion art. sie rg sah ACNET 
of a ; _ 
| » aclT RT ie te ™ . ie 
RTE Lamy et ine" ont nosis pee 
. ‘ is - oii . ’ 
ule gn .ilsay few ia'y sign oD os nolonbovil . 
peter ty ted telAvoEs> ert 2o louite gal New eayrigcse : 
See Cee 
+ ore Ov Iese: To ootaig gies meh? itevy botuisases © ie ie 
wl ‘ 4 7 = 
) ee oreel evidvanet satang kecrtiad Cait bee apebarttot a bras viarle A 
wt on ¢ yale acs ey eel? iw tae at tev _ ei AT ever te 
beac" seth IBRD OF NO eyes 8 to Ye. Ww act by a tHe 
fie: a fe 
igiow one “erttulid Irag" aig ee uM Yieo wee ae drat wera 


aa rage 7 


_ ’ aunnnnscns vin a 


arte A aeoiuiertes ve 


a ‘guahew act oF act et socal re me on — 
a \ i: Pa ed 

AD ager iaes . J ATH, a 4 ool omnia 
ay | crags © & 


L “ gistemmteg | > bewt aid tl otignagey tt vvczieven 
: at 6 vere?) ai aia saruiy 
; Jeet on } = 


This study will focus on a number of areas of shaly sand analysis. Parameters from 
geophysical well logs can be used to discern lithology through the construction of an M-N 
plot. Various porosity logs and the corrected bulk density log are used to develop the 
parameters: volume of shale, volume of dispersed shale, effective porosity and bulk 
volume clay fraction. These parameters are then compared through the use of cross-plots 
with the core-derived parameters to establish relationships between them. 

The resultant relationships between core analysis and the geophysical well logs are 
limited in their application to the Clearwater Formation within a limited geographic area 
(Wyllie and Rose, 1950). Reservoir conditions cannot be assumed to be uniform 
throughout the Cold Lake oil sands deposit. However, the conclusions of the present 


study can serve as a guide for further study of this deposit. 


Function of Well Logging Tools 

It is important to understand the functioning and limitations of the various 
geophysical well logging tools which were run on well 10-—17—-66-5W4 (Figure 28 and 
Appendix |) before analyzing the shaiy sand of the Clearwater Formation (Fertl, 1979). 

The resolution of most geophysical well logging tools is in the order of 0.6 m (2 ft) 
Consequently, some of the more finely interbedded and the interlaminated lithologies will 


not be detected as discrete units by the geophysical well logs. 


Caliper 

The caliper records the diameter of the borehole and indicates areas of 
borehole collapse. The density and neutron readings depend upon a good contact 
between the tool and the borehole wall for accurate readings. The caliper readings 
for well 10-—17-66-5W4 indicate that little borehole collapse has occurred and that 


readings can be considered reliable. 


Spontaneous Potential / Dual Induction Laterolog 

This is an electric log which imparts an electric current into the formation. 
Three simultaneous resistivity recordings were run in well 10-17-66-5W4. The 
laterolog 8 records a shallow reading of the resistivity of the zone invaded by 


drilling fluids. The medium induction log explores further into the formation and may 
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penetrate the undisturbed zone. The deep induction log penetrates furthest into the 
formation and usually records true resistivity. The conductivity is simply the inverse 
of the true resistivity. 

The spontaneous potential log measures the naturally occurring potential 
differences between a surface electrode and a moveable electrode in the mud 
coiumn. The spontaneous potential detects the electrical potential created by the 
difference in the salinities of the formation water and the drilling mud filtrate. The 
presence of an effective shale is indicated by a deflection to the right on the 
Spontaneous Potential curve. The spontaneous potential is most useful as a shaliness 
indicator in water-bearing sands of low to moderate resistivity where the shale is 
laminated. When hydrocarbons and dispersed shales are present, the spontaneous 
potential response is subdued and the caluclated volume of shale values are too high 
(Helander, 1980). Because the Clearwater Formation is saturated with heavy oil and 
consists of dispersed, laminated and structural shales, shaly sand analysis using the 


spontaneous potential log was not performed. 


Gamma Ray Log 

The gamma ray log is considered to be a lithology log which can indicate the 
shaliness of a formation by measuring the intensity of naturally-occurring gamma 
rays. The gamma ray logging tool accompanying the borehole compensated sonic 
log is situated in the centre of the borehole and may be less sensitive due to gamma 
ray attenuation in drilling muds. Consequently, the readings from this log are lower, 
and possibly less reliable, than those from the compensated neutron-formation 
density log. Therefore, the gamma ray values obtained from the compensated 
neutron-formation density log will be considered more reliable and representative. 

Shales are generally more radioactive than sands or carbonates. However, 
the gamma ray log can be misleading in shaly sands since not all clays are radioactive 
(e.g., kaolinite and chlorite) and consequently will not be detected by the gamma ray 
log. Furthermore, the gamma ray log will detect radioactive minerals other than 
those in shales. For these reasons it cannot be considered a completely reliable 


lithologic tool. 


aineres? oe} riruste: ic ; at 

burn estt Wi aborrogle endevcee 

att yo hetaerc ‘oimarog tsahtoale © 

BAT ore th? Gum pri art bans ike 

wtf to TAghy art a oe 

verwwete & cm Wise 720m 2) ier BIS BHreag lad 
uf sterile art Bart yfivitniaet grayabon a5 

“ paaknamermcres we tt Inaaerg ere salete bow seqellp tem age 

ip Oe ay eaualny olerie 19 eenuioy bare auts ae we! ” 
one o-yvanel ide bate utnn ai none | | 
si) Brey. ako yigns Onek yleig: gamete lew toy i 


ort areninn vali 9k yaoi s wc ot peveblanda gal ye erin OFT : ia, 
erieneg geirmpade ay larnan 1 aPorotey ent oni use rapa) a sar _ 
aheere batsaragnan oisdaned wd? CAN ARANOSoS 20 enti col gen kenrg, att aye io 
ghia sa a ction oor sot an pins tl aca ah 
Ae oe iyiet aert rear cerita ti yibraup ean 3 od Opeth 1 siteurr 
aensonrety Porep HL on badmamagagnron ert? ro oot } etdaitey eae | 
aarenireqne: at mont beniatdo eouley yet oe ie 
sitio ts re a ae soho seris “no UUs 
mega sta iaaherdl'y neon ot a ie ning» ck: ‘ BY 
— aie aia 2 weir | ‘ barennrers sah taney 
onan verte ae wrome 
| ‘diteder vie Narn; 


| 5 
a ea 


wed 


Corrected Bulk Density Log 

This tool uses gamma rays to measure the electron density of the formation. 
Its measurements are a function of the density of the rock matrix and the density of 
the fluids and/or gases filling the pore space. It detects both pore water and 
clay-bound water (Juhasz, 1979). The density porosity log converts the bulk density 
into apparent sandstone porosity units, assuming a clean sand density of 2.65 


gms/cc. 


Sonic or Acoustic Log 

The sonic log measures the travel time of sound through the formation. The 
values generated are inversely proportional to the speed of sound and are a function 
of lithology, porosity and the types of fluids present. The sonic porosity is usually a 
good approximation of the total porosity, that is, the effective plus the 
non-effective porosity (Tixier et al, 1959). 

The sands of the Clearwater Formation are not compacted and the sonic log 
signal is decelerated because the grains are not always in direct contact with one 
another. A compaction correction calculated from the sonic log must be applied to 
the sonic porosity. The compaction correction reduces overly optimistic sonic 


porosity values in uncompacted sands. 


Neutron Log 

The neutron log is based on the concept that the mass of a neutron is 
approximately equal to the mass of a hydrogen atom. This log measures the amount 
of hydrogen in the formation by bombarding the formation with neutrons which then 
collide with the hydrogen nuclei (Fertl, 1972). Because it detects the hydrogen 
atoms in both water and in the lattice of the clay minerals (Hill et al, 1979; Juhasz, 
1979), it can be used in the determination of both the volume of shale and the 
porosity of the formation (Truman et al, 1972). It is used extensively in shaly sand 
analysis because it gives good volume of shale values in hydrocarbon-saturated 


sediments (Poupon and Gaymard, 1970). 
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The Effects of Clays on Geophysical Well Log Response 

Shaly sand analysis attempts to deal with the distribution of various types of clays 
and their effects on the geophysical well log reponse. 

The manner in which clay is distributed in the sediment affects the geophysical well 
log response. Three forms of clay distribution are recognized (Fertl, 1972; Heslop, 1975; 
Ransom, 1977; Helander, 1980). They are dispersed clay, structural clay and laminar clay. 
The core and thin section studies show that all three forms are present in the Clearwater 
Formation. 

Dispersed clay can be detrital or authigenic and is found on grain surfaces or 
occupies the void space between grains. This habit of shaie affects both porosity and 
permeability. The shale acts as a permeability barrier for the transmission of fluids. The 
framework porosity is decreased by the combined effects of the volume of dispersed 
clay and the bound or adsorbed water associated with the clay. 

Structural clay supports overburden pressure. It is composed of fragments or 
crystals which are part of the rock framework. While this type of clay does not affect the 
total pore volume, it does reduce the effective porosity by the amount of adsorbed water 
associated with the clay. 

Laminar clay consists of both dispersed and structural clays. It occupies space in 
both the rock framework and in the void space between grains. Since it contains both 


dispersed and structural clays, it reduces both porosity and permeability. 


Clay Composition 

Each clay type affects the geophysical well log response differently. The 
water contained in the pores of a shaly sand is closely associated with the clay as 
adsorbed water. The chemical composition of the clay and the amount of adsorbed 
water associated with it causes deviations from the geophysical response as 


compared to that of a quartz sand. 


Kaolinite 
Kaolinite is non-radioactive and is detected as a sand by the gamma 
ray, resistivity and spontaneous potential logs (Hilchie, 1982). It is therefore 


termed a non-effective clay. The gamma ray and spontaneous potential 
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responses for this clay are low and the resistivity response is high. The 


density of kaolinite is 2.43 gms/cc which is lower than that of quartz sand. 


Illite 

Illite is a radioactive effective clay and is detected as a shale. Both the 
spontaneous potential and gamma ray logs record illite as a shale. The density 
of illite is very close to that of quartz sand and is 2.61 gms/cc. 

The presence of this clay is indicated when the neutron porosity is 
much higher than the density porosity. The density porosity is lower than the 
true porosity when illite is the major constituent in the sediment. Iilite cannot 


be differentiated from smectite on the sonic porosity log. 


Smectite 

Smectite is usually non-radioactive, but can be a radioactive effective 
clay when uranium in solution is attracted to the charged smectite clay 
particles. When uranium is present in large amounts, smectite is detected by 
a rapid increase in the gamma ray log. The presence of smectite is also 
indicated by low resistivity values. The density of this clay is 2.30 gms/cc 
which is lower than the density of quartz sand. Consequently, sediments that 
contain large amounts of smectite will register low bulk density values. 

All three porosity logs will record much higher porosities than the 
true porosity when smectite is present. While all clays exhibit interparticle 
swelling, smectite exhibits both interparticle and interlayer swelling and 
provides the largest total surface area for adsorption of water (Jenkins and 
Bush, 1971). Smectite can have up to four layers of adsorbed water in 
uncompacted sands, such as those of the Clearwater Formation. The neutron 
porosity log will be affected by the number of layers of adsorbed water 
since it detects hydrogen atoms (Hill et al, 1979). While all porosity log 
values will be much higher than the true porosity, the neutron log will exhibit 


much higher porosity values than the others when smectite is present. 
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Chlorite 

Like kaolinite, chlorite is a non-radioactive non-effective clay which is 
detected as a sand by the gamma ray and spontaneous potential logs. It has a 
low gamma ray and high resistivity response and a density of from 2.78 
gms/cc to 3 gms/cc which makes it slightly more dense than a quartz sand. 
This high density is reflected by the density porosity log which is usually very 
low as compared to the sonic and neutron porosity logs. The sonic and 
neutron porosity logs exhibit higher porosities similar to those for kaolinite. 
Consequently it is difficult to differentiate between kaolinite and chiorite on 


geophysical well logs. 


Summary 

X-ray diffraction analyses of samples from well 10-—17-66—-5W44 indicate that 
kaolinite, illite, smectite and chlorite are present in varying amounts with the sediments of 
the Clearwater Formation. These clays are found in dispersed, structural and laminar 
forms. The mixture of effective and non-effective clays (Hilchieé, 1982) and the wide 
variety of forms these clays take (Ransom, 1977), complicates the shaly sand analysis to 
the extent that clear linear relationships between geophysical parameters and core-derived 


parameters cannot be expected. 


C. Shaly Sand Plots 

The data for the various plots are derived from the hand-digitization of the 
geophysical well logs, the mathematical transformations of these data and the core 
analysis values. In the case of the M-N plot, all 301 points digitized were used. 

in cross-plots using the data derived from the geophysical well logs, the points 
representing the calcareous cemented zones were omitted because they obscure the shaly 
sand analysis. The number of points plotted is reduced to 264. 

The cross-plots that incorporate data derived from the core analysis do not include 
thick sections of shale (all of stratigraphic units A and J and the upper and basal shales of 
Stratigraphic unit D). in these cases, core was not sampled and the number of points 


plotted is 165. 
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Sources of Error in the Determination of Core-derived Parameters 

Sampling procedures and the method of obtaining weight percentages of oil and 
water contribute greatly to the sources of error in the core analysis of shaly sands. 

A considerable amount of water can be added to the samples during coring so that 
weight percent water becomes higher than that found in situ. Core porosities as high as 40 
percent are probably the result of water gain through core expansion and imbibition during 
coring (Eade, 1975). Some of the core porosities determined for the Clearwater 
Formation are as high as 43 percent. Conversely, water loss due to evaporation after 
coring can range up to 60 percent (Eade, 1975). This affects the values for the weight 
percent water and the core porosity, although the weight percent oil is not affected. 

Higher core porosities can also be caused by core expansion due to the release of 
overburden pressure (Woodhouse, 1976; Collins, 1977). The amount of expansion has 
been estimated to be 6 percent in the Athabasca oil sands area (Woodhouse, 1976). While 
no figures are available for the Clearwater Formation in Cold Lake, the amount of 
expansion is probably higher due to the greater thickness of overburden tn this area. Core 
expansion further complicates the relationship between core- and log-derived parameters. 
The Dean-Stark method (Kidston, personal communication) was used to analyze the oil 
sands samples. In this method, toluene is used to separate the oil from the sample. 
Because toluene boils at Ill°C., both adsorbed water and moveable water are removed 
from the sample. This contributes to overly-optimistic values for core porosity (Bush and 
Jenkins, 1969; Jenkins and Bush, 1971; Woodhouse, 1976}. The Dean-Stark method 
assumes that fiuids occupy the total pore space and that there is no air present in the 
sample. This contributes to inaccurate values for weight percent bitumen and water. 

These sources of error make it difficult to quantitatively resolve core-derived with 
log-derived parameters, although general trends can be observed. It is imperative that 
better methods of core extraction and more precise methods of core analysis be 


developed to minimize these problems. 
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M-N Plot 

The M-N plot (Figure 29) is used te determine general lithology using the neutron, 
density and sonic logs (Pirson, 198 1b). The values of M and N (Appendix Ill) are dependant 
on fluid and matrix characteristics, but are independant of porosity (Schlumberger, 1979). 
The plot is based on the concept that the various well logging tools respond in different 
ways to the various material components of a complex lithology (Helander, 1980). The 
logging devices have a characteristic response to the major components of reservoir 
rocks: quartz, calcite, dolomite, anhydrite and gypsum. 

When the resultant plot is compared with the M-N plot for points of single mineral 
formations, the shaliness of the Clearwater Formation is apparent. All points piot well 
below the clean sand point and are shifted in the direction of the shales. 

The sample points were characterized as calcareous-cemented sands, relatively 
clean sands, shaly sands or shales by the dominant lithology found at equivalent depths in 
the core. When these categories are plotted on the M-N plot the shaliness of the 
formation can be observed influencing all lithologies by pulling them down and to the left. 
There is some separation of shaly sands and shales into two end members, but there is 
also a large area of overlap. The shale category represents the shales in stratigraphic units 


A and J. 


Corrected Sonic Porosity versus Neutron Porosity 

The sonic porosity was derived mathematically (Appendix Ill). A compaction 
correction was applied to the sonic porosity because the overburden pressure is 
considerably less than the formation pressure (Tixier et al, 1968). The lack of compaction 
of the sands causes considerable attentuation of the sonic signal. The need for a 
compaction correction is indicated when the sonic travel time in shales closely associated 
with the sands exceeds 100 microsec./ foot (Tixier et al, 1959). In the case of the 
Clearwater Formation in well 10-—17-66—-5W4, the closely associated shales of 
stratigraphic unit D exhibit a travel time of approximately 130 microsec/ foot. 
Consequently, compaction correction of 1.3 was applied to the sonic porosity. 

All of the log-derived porosities plot in a narrow band between 30 percent and 60 


percent. Porosities are high due to the effects of the clays. The corrected sonic 
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a'e region 


M versus N lithoplot, showing the effect of 


clays on the lithology of Clearwater sediments. 
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porosity-neutron porosity plot (Figure 30) exhibits the best direct relationship (r=0.63) 
since both logs reflect total porosity. However, the neutron log also detects hydrogen in 
the clay crystal lattice. This would account for the average 6.5 porosity unit increase in the 
neutron porosity as compared to the sonic porosity. The discrepancy between the neutron 
and sonic logs is particularly evident in the lower porosities where the sample points are 
more shaiy. In these areas the neutron porosity can be up to 19 porosity units higher than 


the sonic porosity. 


Shaliness Indicators 

Several plots were generated to obtain more realistic values for volume of shale 
and effective porosity throughout the Clearwater Formation. Shaliness indicators in the 
porosity logs, the corrected bulk density log and the gamma ray log give good values or 
the upper limit for the volume of shale in a shaly sand formation where there are no clean 
100 percent water-saturated sands which can serve as areference. Usually the lowest of 
the upper limits for volume of shale values are taken as the most realistic values (Poupon 
and Gaymard, 1970; Helander, 1980). 

Effective porosity calculated from multiple log analyses provides the most 


meaningful values for true porosity of shaly sands (Heslop, 1975). 


Corrected Density versus Neutron Porosity 

This cross-plot was generated to obtain values for the volume of shale and 
effective porosity. This method was first described by Poupon and Gaymard (1970) 
and is discussed by others (Krug and Cox, 1976; Ransom, 1977; Helander, 1980). A 
mathematical approach to the determination of volume of shale and effective 
porosity was described by Krug and Cox (1976). 

In the corrected bulk density-neutron porosity treatment (Figure 31) the wet 
clay point is defined in an associated shale with the highest neutron porosity. This 
assumes that the wet clay represents a shale which is 100 percent effective. If Is is 
not 100 percent effective, volume of shale figures will be high for the more shaly 
points, but close in the cleaner sands (Hilchie, 1982). The clean quartz sand point is 
placed at a neutron sandstone porosity of 0 percent and a density of 2.65 gms/cc 


which is the density of a pure quartz sand. The water point is placed at 100 percent 
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Corrected bulk density versus neutron porosity, 
illustrating the construction of the shaly sand 
plot. 
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Detail of shaly sand plot, showing relationship 
of points to volume of shale and effective porosity. 
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porosity and a water density of 10 gm/cc. This trianguiar plot defines end points 
for a pure quartz sand, a pure shale and 100 percent porosity. The points on the 
corrected bulk density-neutron porosity cross-plot are projected onto these lines 
and values for volume of shale and effective porosity are obtained (Figure 32). 
Points which plot outside the triangle of the shaly sand plot represent anomolous 
shales which have values for bulk density and neutron porosity outside the range of 
accepted values. 

Several factors influence the shaly sand plot generated by the corrected bulk 
density-neutron porosity cross-plot. The influence of the density of large amounts 
of heterogeneous solid matter can shift the plot (Figure 33). Several components of 
the sediment in the Clearwater Formation have an effect on the density log. As can 
be seen, values for the corrected bulk density of the digitized points range between 
20 gms/cc and 2.40 gms/cc which is much lower than that of sand which is 2.65 
gms/cc. Biotite, zircon, siderite, hematite, magnetite and pyrite are all much more 
dense than a pure quartz sand and tend to pull the points down (Edmundson and 
Raymer, 1979). Given the low densities of the digitized points, it seems that the 
more dense minerals are not present in sufficient amounts to significantly affect the 
bulk density. With the exception of chiorite which is slightly more dense that quartz, 
the clays identified and the carbonaceous debris lower the density values. These 
components appear to have the most effect on the density log in the formation 
under study since the bulk density values are low when compared with the density 
of aclean quartz sand. 

The second factor to take into account is the effect of light hydrocarbons. 
Corrections must be made for the presence of light hydrocarbons in the formation 
since they affect the readings of the sonic, density and neutron logs (Poupon and 
Gaymard, 1970). The degree of influence of the hydrocarbons is a function of the 
degree of hydrocarbon saturation and the density of the hydrocarbons (Juhasz, 
1979) and is determined by the calculation of the hydrocarbon index. The neutron 
porosity of the hydrocarbon can then be caiculated and if it is significantly different 
from 100 percent a hydrocarbon correction is indicated. The density of the 


hydrocarbons present in the Clearwater Formation was determined by several 
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independant studies to be 0.999 gms/cc at I5°C. which is very close to the density 
of water. The neutron porosity was determined to be 100 percent for the 
hydrocarbon. This means that the hydrocarbons have no effect on the neutron, 
sonic and density logs. 

The presence of light hydrocarbons was also checked for by the 
neutron-density overlay method (Tixier et al, 1968) which showed that no gas was 
present in the formation. 

Consequently, no correction was needed for the vaiues of volume of shale 
and effective porosity. 

The corrected bulk density-neutron porosity cross-plot shows that no clean 
water saturated sands are present in the formation. The concentration of points at 
the clay point indicates that the shales in the formation have a high clay component 
and a low silt component (Bosworth, personal communication). The volume of shale 
ranges from 25 percent in the cleanest sands to 100 percent in the pure shales. 
Effective porosity ranges from 0 percent in the shales to 31 percent in the cleanest 
sands. This cross-plot gave the lowest upper limit for the volume of shale and was 
used in further computations. 

A cross-plot was generated to compare the effective porosity and the 
volume of shale (Figure 34). It exhibits a good inverse relationship (r=-0.9 1), 
showing that as the volume of shale increases the effective porosity decreases. This 
relationship was also noted by Heslop (1975). 

The effective porosity also exhibits a good direct relationship (r=0.91) with 
the density porosity (Figure 35). Although there is some scatter, in general the 
porosity is much higher than the effective porosity at lower porosity values 
(representing shales) and converges at approximately 50 porosity units in the cleaner 


sands. 


Apparent Sandstone Density Porosity versus Neutron Porosity 
This cross-plot was generated to compare and confirm the values for 
volume of shale and effective porosity obtained by the cross-plot of corrected bulk 


density and neutron porosity. 
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Because the density and neutron tools detect both effective and 
non-effective porosity, amethod was developed to calculate effective porosity 
only. This shaly sand plot (Figure 36) was first described by Poupon et al (1970) and 
is discussed in other papers (e€.g., Ransom, 1977; Helander, 1980; Hilchie, 1982). 

In the density-neutron porosity treatment, a wet clay point is defined in an 
associated shale as 100 percent volume of shale and 0 percent effective porosity. 
The position of the wet clay point is dependant upon the type of clay, the amount of 
water associated with the clay and the mineralogic response of the neutron tool 
(Ransom, 1977). In this case, the clay point is defined by a smectite-dominated shale 
because it has high density and neutron porosity values. This is confirmed by x-ray 
diffraction. Points which plot outside the triangle of the shaly sand plot represent 
anomolous shales which have values for density porosity and neutron porosity 
cutside the range of accepted values. A clean water sand line is constructed from 
the O percent neutron and density porosities to the fluid point at 100 percent 
neutron and density porosities. These two lines are then divided from 0 to 100 
percent and represent volume of shale and effective porosity values. Plotted points 
are projected onto these lines and values for volume of shale and effective porosity 
are obtained (Figure 37). 

The density-neutron porosity cross-plot confirms that there are no clean 
water-saturated sands in the formation. Volumes of shale range from 20 percent in 
the cleaner sands to 100 percent in the shales. Effective porosity ranges from 0 
percent in the shale to 30 percent in the cleaner sands. The values for volume of 
shale and effective porosity compare well with those generated by the corrected 
bulk density-neutron porosity cross-plot but are slightly less optimistic. 


Consequently, these values were not used for further computations. 


Gamma Ray Index versus Volume of Shale 

The gamma ray log from the compensated neutron-formation density was 
used to obtain gamma ray index values (Appendix Ill). The gamma ray index values are 
an indication of the volume of shaie (Woodhouse, 1976} and can be compared to 
the volume of shale generated by the corrected bulk density-neutron porosity 


cross-plot (Figure 38). 
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FIGURE 36. Density porosity versus neutron porosity, 
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illustrating the construction of the shaly 
sand plot. 
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FIGURE 37. Detail of shaly sand plot, showing relationship 


of pointsto volume of shale and effective porosity. 
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GURE 38. Gamma ray index versus volume of shale, showing 


consistently larger percentage of shale indicated 
by gamma ray index. (r = 0.85) 


oo ¢ 


oo ° 


0° bd 


Volume of Shale (%) 


FIGURE 39. Bulk volume clay fraction versus volume of shale, 


illustrating that the bulk volume of clay and the 
adsorbed water on the clays are directly related 
to the volume of shale. (r = 0.82) 
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This cross-plot exhibits a good indirect relationship (r=0.85) between these 
two values for the volume of shale. The gamma ray index values are consistently 
higher than those from the corrected bulk density-neutron porosity cross-plot. This 
suggests that the gamma ray log is not simply detecting the shale radioactivity, but 
also radioactive minerals in the sands. The presence of glauconite and feidspars in 
the sediments can greatly increase the gamma ray readings (Edmundson and Raymer, 
1979). This will create higher apparent volume of shale values. Consequently, the 
gamma ray index-generated values for volume of shale were not used in further 


computations. 


Relationships between Log- and Core-derived Parameters 

Both the geophysical log and core data were used to generate a number of 
parameters. These values were then compared. 

The high variablity of clay composition and form and the previously-discussed 
problems associated with core analysis meant that few relationships between these 
parameters could be discerned. The following is a presentation of the clearest 


relationships found for the Clearwater Formation in well 10-—17-66-5W4. 


Bulk Volume Clay Fraction versus Volume of Shale 

This cross-plot (Figure 39) attempts to resolve core-derived and 
log-derived parameters. The bulk volume clay fraction represents the 
fraction of the bulk volume occupied by adsorbed water on the clays and is 
calculated by subtracting effective porosity values from the core porosities 
(Appendix Ill). The plot indicates that the bulk volume of clay and 
consequently the percentage of adsorbed water is directly related (r=0.82) 
to the volume of shale, especially in the cleaner sands. This confirms that the 


volume of shale is a gocd total shaliness indicator. 


Volume of Dispersed Shale versus Volume of Shale 
The volume of dispersed shale in the pore space is calculated by 
sonic and neutron porosity measurements (Appendix III). The cross-plot 


(Figure 40) shows that as the volume of shale increases the total volume of 
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FIGURE 40. Volume of dispersed shale versus volume of shale, 
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Calculated Effective Porosity (%) 
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illustrating a direct relationship and the 
presence of laminar and structural shale. (r = 0.79) 


Effective Porosity (%) 


FIGURE 41. Calculated effective porosity versus effective 


porosity, illustrating the relationship between 
log-derived and core-derived porosities. (r = 0.75) 


-atedz to smlov 
ai? bes qin 
(RTO 9) atone Toes. ns 


oF 


dispersed shale also increases (r=0.79). However, the volume of shale is up 
to twice the amount of dispersed shale. This confirms that considerable 


amounts of structural and laminar shales are present in the sands. 


Calculated Effective Porosity versus Effective Porosity 

This cross-plot attempts to resolve all log-derived porosity 
measurements and the core-derived porosity measurement. The calculated 
effective porosity is a function of the corrected sonic porosity, the density 
porosity and the core porosity (Appendix Ill). The effective porosity is a 
function of the corrected bulk density and the neutron porosity. 

The calculated effective porosity is consistently higher than the 
effective porosity (Figure 41), especially in the lower porosity values 
(r=0.75). The two independantly calculated porosities tend to converge at 
higher porosity levels. The major reason for the discrepancy is that the core 
porosities are overly inflated because the clay porosity measurement 
reflects adsorbed water as well as free pore water. The core porosity 
measurements may also be infiated due to core expansion and other 


previously-discussed sources of error. 


Volume of Dispersed Shale versus Effective Porosity 

This plot (Figure 42) shows that two independant computations of the 
three log-derived porosity measurements give correlatable results. This 
reflects the internal consistency of the log-derived parameters. The 
cross-plot exhibits a good inverse relationship (r=-0.84), showing that as the 
volume of dispersed shale increases, the effective porosity decreases. 

The plot can be divided on the basis of lithology. There is a good 
concentration of points in the shaly sands which contain between 12 percent 
and 25 percent dispersed shale and between 20 percent and 30 percent 


effective porosity. 


Cae 


yneceg Paves eee 
befelioliss ont trace weaen: is anc i 
Planad ant .yiieo1wg DHS berate wt 


si yreow, ovisoet te ort Mb 


gr? merit senpey (inetaianen ei e bi 
gouwtey VNeotog weil diab eben: - ' "4 ie 
i By eves -d nas soitecnngeattlaniet eo 
Svs arts fer't 2) oneqe nals oat WF dose 10 > 
tre ne wegen yOu Yay ae we phn: te | 
yilaarod e105 rit . iatew sod eet eran not aeannnetanesen 
verte brig oaienn qs ered of seb oatnitnl ad oaln yom sinemawanain | 
Lone bo senwEs beseloeb: vienna 


va 
eancrefabe 9 tua aa ane a sel Pe 

apt To-noctinn yee: 7 tan riod: chun tert aworte (GA an 
(aT alu SeiaaBIE NCD Vig eremenuBoN yieOTS sow 
ent. craton \ty bewitayrgel eit 10 (retalane tae f toatOn 

arid ony ser grtiwerts if. O- “= yéjilancioniey seve Reo aiid lense, | 
dvesprosh wieo req ovitaette ert aseso70ri Gis guremonk 
boon ei ered ygotorttil te aifed art! no bativil ea palo is 


=m T° 
ch nl 


mre i maw inbnnaiameaiens Yelwareatt i i nti mgt 


fn 


ra) ae 
io wea 3 


Porosity  (%) 


Core 


(decimal %) 


Effective Porosity (%) 


FIGURE 42. Volume of dispersed shale versus effective 
porosity, showing decrease in porosity with 
increase in dispersed shale. (r = -0.84) 


Density Porosity  (%) 


FIGURE 43. Core porosity versus density porosity, 
illustrating the poor relationship between 


log-derived and core-derived porosities. (r = 0.55) 
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Core Porosity versus Apparent Sandstone Density Porosity 

All log-derived porosities were cross-plotted with the core porosity. 

The relationship between the core porosity and the apparent 
sandstone density porosity (Figure 43) does not show a very consistent 
relationship (r=O.55) and much scatter is evident. The difference between 
the core and density porosities is greatest at lower porosity levels, but 
values tend to converge at higher levels. The discrepancies between these 
two porosity measurements are caused by a combination of highly variable 
mineralogic control on the density porosity and the inflated core porosity 


measurements. 


Core Porosity versus Corrected Sonic Porosity 

The sonic porosity is considered to be the best approximation of 
total porosity (Juhasz, 1979). Total porosity is also represented by core 
porosity. The cross-plot (Figure 44) indicates a fair direct relationship 
(r=0.66) between corrected sonic porosity and core porosity. On average, 
the corrected sonic porosity is 3.2 porosity units higher than the core 
porosity, although it ranges between 5 porosity units lower than the core 
porosity at one point to 7 porosity units higher at several points. These 
discrepancies are primarily due to the fact that corrected sonic porosity was 
measured every 0.15 m (6 in), whereas core porosities were generalized 
from one plug sample to represent core up to! m long. The piot suggest that 
a compaction correction of between 1.3 and 1.4 will reconcile the core and 


sonic porosities in Clearwater Formation shaly sands. 


Core Porosity versus Effective Porosity 

One major aim of shaly sand analysis is to link effective porosity and 
core porosity in a quantitative fashion. Since total core porosity includes 
adsorbed water on the clays and effective porosity does not, the difference 
between the two values should indicate the amount of bound water (Heslop, 


1975; Ransom, 1977; Juhasz, 1979). 
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FIGURE 44. Core porosity versus corrected sonic porosity, 
illustrating a good relationship between log- 
derived and core-derived porosities. (r = 0.66) 
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FIGURE 45. Core porosity versus effective porosity, 
illustrating the large discrepancy between 
these parameters. (r = 0.59) 
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The plot of core porosity versus effective porosity (Figure 45) 
indicates that all samples contain adsorbed water (r=0.59). The core 
porosities are consistently much higher than the effective porosities. The 
proportion of adsorbed water is greater in the region of lower effective 
porosity values which represent the more shaly sands. The effective and 
core porosities are less divergent in the cieaner sands which contain less 
clay. This is indicated by the convergence of the line of equal porosity values 
with the line defining the reiationship between the core and effective 
porosities. 

The cleaner sands plot in the upper right-hand part of the cross-plot. 
Here the effective porosity ranges between 7 porosity units and 17 porosity 
units lower than the total or core porosity. This discrepancy is caused by a 
combination of greater amounts of adsorbed water in the clays of the more 
shaly sands and the previously-discussed inaccuracies in obtaining core 
porosities. On the average, the effective porosity in the cleaner shaly sands 


is 12 porosity units lower than the core-derived porosity. 


Effective Porosity versus Bulk Volume Water 

The plot of the effective porosity versus the bulk volume water 
(Figure 46) illustrates the problems with the core analysis techniques 
discussed previously. The wide degree of scatter (r=-0.24) is a function of 
the generalized values obtained by the core analysis. 

The general trend is to increased bulk volume water in the shaly sands 
with lower effective porosity values. This is due to the greater amounts of 


adsorbed water in the more shaly samples. 
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FIGURE 46. Effective porosity versus bulk volume water, 
showing a poor relationship and reflecting an 
increased amount of bound water in shalier 
sands. (r = -0.24) 
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Vi. Conclusions 


A. Sedimentologic implications 


ie 


The facies association analysis and grain size analysis indicate that the Clearwater 
Formation consists of two coarsening-upward sequences. The first represents the 
deposition of poorly developed offshore marine bars and ridges on a muddy shelf 
base. The second represents the deposition of nearshore sediments, including tidal 
bundles, storm graded deposits small-scale bars and associated interbar muds and 
silts. This second sequence also contains a shale deposit which represents a period 
of non-deposition of coarser-grained sediments. 

This environment of deposition is highly complex and better well control is 
needed in order to trace the position of the sand bars and ridges. The site holds 
promise for the discovery of more potentially hydrocarbon-bearing fine-grained 
sand bars and ridges. These bars will be of limited vertical and horizontal extent. The 
silty shale and shale interlaminae and interbeds may at times act as permeability 
barriers and serve to trap hydrocarbons. However, they also may complicate the 
recovery process. Many of these laminae and beds are discontinuous, rather than 
continuous, which will be advantageous to recovery. 

Calcareous cemented zones may also act as permeability barriers, depending 
upon their aerial extent. 

The Clearwater Formation in the study region consists of very fine-grained 
feldspathic litharenites and litharenites. The mineralogy of this formation is complex 
and highly variable. 

In situ hydrocarbon recovery methods will have to be deveioped that can 
deal with the large amounts of complex authigenic clay mineral assemblages in order 
to ensure that the reservoir will not be destroyed by clay alteration and migration 


into the pore spaces. 
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B. Geophysical Well Logging Implications 


i 


The corrected sonic porosity is the best indicator of total porosity (effective and 
non-effective) and correlates best with the core porosity. On the average, the 
corrected sonic porosity is 3.2 porosity units higher than the core porosity. A 
compaction correction of the sonic porosity of between |.3 and |.4 is indicated in 
order to reconcile sonic and core porosities. 

The cross-plot of corrected bulk density and neutron porosity gives the most 
realistic effective porosity values and is a good shaiiness indicator. The cleanest 
sand in the Clearwater Formation is found in Facies F and has a volume of shaie in the 
order of 25 percent. The highest effective porosity is 31 percent. The high 
percentage of shale indicates that the formation is at the borderline of producibility. 
The actual effective porosity of the cleaner shaly sands is approximately 12 
porosity units lower than that indicated by core analysis of well 10-17-66-5W4. 
The proportion of dispersed shale increases in the more shaly sands. The dispersed 
shale generally accounts for one-half of the total shale present. The remainder is in 
the form of structural and laminar shale. - 

The gamma ray is not a good indicator of shaliness in shaly sands which contain large 
amounts of non-effective clay, such as kaolinite. It is also not a good shaliness 
indicator in sediments which contain radioactive minerals, such as glauconite and 
certain feldspars which cause high gamma ray values, not reflective of shale. 

Better shaly sand analysis techniques are needed to deal with formations, such as the 
Clearwater Formation, which contain both effective and non-effective clays and a 
mixture of dispersed, structural and laminar clays. A natural gamma ray 
spectroscopy log would help differentiate between clay / shale types as would a 
density lithology log. These logs may also be abie to identify other radioactive 
minerals such as glauconite. Other computer programs developed by the well 
logging industry may be able to establish approximate sand, silt, shale and matrix 
percentages. 

Core analysis techniques should be refined to deal more effectively with core in 
unconsolidated shaly sands. In particular, a core analysis technique Is required that 


more clearly distinguishes between adsorbed water on clays and grain surfaces and 
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Plates 
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Plate 1: Core Photographs for Stratigraphic Units J to H 


1. Stratigraphic unit J. 
Facies Association 2. Heavily Bioturbated Silty Shale to Silt and Sand. 
Horizontally laminated; interlaminae and interbeds of lightly oil stained argillaceous 
silt; bioturbated by traces of Planolites nicholson (1873). 
Syndepositional fault at 475 m (1558.25 ft) marks contact with stratigraphic unit |. 
Well 10-17-66-5W4: 475 m (1557.5 ft) at left to 475.8 m (1560 ft) at right. 
Scale: 1 cm = 7.6 cm 


2. Stratigraphic unit |. 
Facies Association 2. Heavily Bioturbated Silty Shale to Silt and Sand. 
Horizontally interlaminated to interbedded; heavily bioturbated by traces of Planolites 
nicholson (1873). 
Thicker sandstone bed at 472 m(1548.5 ft) marks contact with stratigraphic unit H. 
Well 10-17-66-5W4: 471.7 m (1547.5 ft) at top left to 475 m (1557.5 ft) at 
bottom right. 
Scale: 1 cm = 7.3 cm 


3. Stratigraphic unit H. 

eee Association 3. Moderately to Heavily Bioturbated Sand with Silt and Silty 
hale. 

Horizontally interbedded to interlaminated; moderately bioturbated by traces of 
Planolites nicholson (1873). 
Large trace of Skolithos linearis haldeman (1840) from 470.5 m (1543.5 ft! to 
470.7 m (1544.4 ft) 
Well 10-—17-66-5W4: 468.6 m (1537.5 ft) at top left to 471.7 m (1547.5 ft) at 
bottom right. 
Scale: 1 cm = 7.3cm 
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Plate 2: Core Photographs for Stratigraphic Units G to E 


1. Stratigraphic unit G. 
Facies Association 3. Moderately to Heavily Bioturbated Sand with Silt and Silty 
Shale. 
Horizontally bedded; silty shale ripple drapes common; lightly bioturbated by traces 
of Planolites nicholson (1873). 
Silty shale intraclasts at 464.4 m (1523.6 ft). 
Well 10-17-66-5W4: 463.3 m (1520 ft) at top left to 465.7 m (1528 ft) at 
bottom right. 
Scale: 1 cm = 6.8 cm 


2. Stratigraphic unit F. 
Stratigraphic subunit F, from 451.1 m (1480 ft) to 452.3 m (1484 ft). 
Facies Association 6. Sand and Carbonaceous Sand. 
Horizontally laminated very fine-grained sand. 
Stratigraphic subunit F, from 452.3 m (1484 ft) to 453.5 m (1488 ft). 
Facies Association 4. Sand, Silt and Silty Shale. 
Horizontally laminated and ripple cross-laminated sands; horizontally and undulating 
silty shale; silty shale ripple drapes. 
Rare bioturbation by traces of Planolites nicholson (1873) 
Well 10-17-66-5W4: 451.1 m (1480 ft) at top left to 453.5 m (1488 ft) at 
bottom right. 
Scale: 1 cm = 8.5 cm 


Stratigraphic unit E. 

Facies Association 4. Sand Silt and Silty Shale. 

Finely interlaminated to very finely interlaminated (note: curved edges caused by 
coring). 

Thick sand at 464.5 m (1524 ft) marks contact with stratigraphic unit F. 

Shale at 462.4 m (1517 ft) marks contact with stratigraphic unit D. 

Well 6-5-66-5W4: 462.2 m(1516 ft) at top left to 465.2 m (1526 ft) at bottom 
right. 

Scale: 1 cm = 8.5 cm 
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Plate 3: Core Photographs for Stratigraphic Units D to B 


1. Stratigraphic unit D. 
Stratigraphic subunits D; and D, from 446.2 m (1464 ft) to 446.8 m (1466 ft) and 
from 447.8 m (1469.2 ft) to 448.7 m (1472 ft). 
Facies Association 1. Heavily Bioturbated Silty Shale. 
Fissile silty shale; remnant horizontal lamination. 
Stratigraphic subunit D, from 446.8 m (1466 ft) to 447.8 m (1469.2 ft). 
Facies Association 3. Moderately to Heavily Bioturbated Sand with Silt and Silty 
Shale. 
Very fine-grained sand and silty shale interlaminated to interbedded. 
Entire subunit bioturbated by traces of Planolites nicholson (1873). 
Well 10-17-66-5W4: 446.2 m (1464 ft) at top left to 448.7 m (1472 ft) at 
bottom right. 
Scale: 1 cm=6.8 cm 


2. Stratigraphic unit C. 
Facies Association 5. Graded Sand, Silt and Shale with Bouma Sequences. 
Interbedded sand, silty shale and argillaceous silt; horizontal laminae, ripple 
cross-laminae; silty shale ripple drapes. 
Complete graded-storm turbidity current deposit with Bouma sequence from 443.6 
m (1455.2 ft) to 443.7 m (1455.7 ft). 
Sand at 441.7 m(1449.2 ft) marks contact with stratigraphic unit B. 
Well 10-17-66-5W4: 444.4 m(1448 ft) at top left to 443.8 m (1456 ft) at 
bottom right. 
Scale: 1 cm =6.8 cm 


3. Stratigraphic unit B. 
Facies Association 6. Sand and Carbonaceous Sand. 
Horizontally laminated and inclined laminae; undulating silty shale laminae. 
Shale at 456.7 m (1498.2 ft) marks contact with stratigraphic unit C. 
at 6-3-66-5W4: 454.8 m (1492 ft) at bottom right to 457.8 m (1502 ft) at top 
left. 
Scale: 1 cm= 7.3 cm 
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Plate 4: Mineralogy 


1. Authigenic feldspar overgrowth (appears as one grain under plain light). Dust at 
boundary. Entire grain is subangular and exhibits abrasion, indicating formation during | 
a previous depositional cycle. Grain surrounded by matrix material. 
Sample 80-4-25; Facies D; 445.3 m (1461 ft) 
Magnification: 390.63 X 


2.  Detrital plagioclase feldspar exhibiting crystal zoning and albite twinning. 
Sample 80-9-5; Facies C; 442.5 m(1451.7 ft) 
Magnification: 390.63 X 


3.  Detrital microcline feldspar exhibiting typical grid twinning. Offset fracture caused 
by overburden pressure. 
Sample 80-9-2; Facies B; 440.2 m (1444.3 ft) 
Magnification: 390.63 X 


4. Volcanic rock fragment; untwinned zoned feldspar microphenocryst within 
finely-crystalline volcanic groundmass. 
Sample 80-9-2; Facies B: 440.2 m (1444.3 ft) 
Magnification: 390.63 X 
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Plate 5: Mineralogy 


1. Polycrystalline quartz of metamorphic origin. Note sutured texture of grain. 
Sample 80-9-4; Facies B; 441.8 m (1449.5 ft) 
Magnification: 390.63 X 


2. Chlorite grain (identified by ultra blue birefringence). Grain is embedded in a larger 
glauconite grain exhibiting soft sediment deformation. 
Sample 80-9-5; Facies C; 442.5 m (1451.7 ft) 
Magnification: 390.63 X 


3. Glauconite replacement of original oolite. Grain has been abraded. 
Sample 80-9-6; Facies C; 442.6 m (1452 ft) 
Magnification: 390.63 X 


4.  Calcareous cement: cone-in-cone structure (to the right) displaces the original fabric; 
calcareous cement (to the left) replaces original matrix, pore spaces and some 
detrital grains. 

Sample 75-1-2; Facies F; 454.5 m(1491 ft) 
Magnification: 156.25 X 
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Plate 6: Mineralogy and Textures 


iG caene veins (white) displacing poorly sorted very fine-grained sand with a slight 
offset. 
Organic debris (black) exhibiting preferred orientation aligned parallel to partings. 
Note that organic debris is also found within the calcite veins. Detrital grains include 
subangular quartz with some chert. 
Individual laminae poorly sorted. 
Sample 80-9-23; Facies F; 454.3 m (1490.4 ft) 
Magnification: 156.25 X 


Ze Organic fragment exhibiting original cell structure. Note larger size of fragment as 
compared to other detrital grains. This is due to the lower specific gravity of the 
organic fragment. Quartz grains at 4 o'clock with respect to the organic fragment 
cemented by quartz overgrowths. Poorly sorted nature of Clearwater Formation is 
also exhibited. 

Sample 80-9-4; Facies B; 441.8 m (1449.5 ft) 
Magnification: 156.25 X 
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Plate 7: Mineralogy 


1. Muscovite mica grain (just above centre of photograph) and sedimentary rock 
fragment (to left of mica) deformed by compaction between quartz grain on top and 
subangular chert grains on bottom. Note rims of oil-stained clays lining the pore 
spaces. 

Sample 80-9-2; Facies B; 440.2 m (1444.3 ft) 
Magnification: 156.25 X 


2. Large mass of authigenic pyrite. Many cubic crystals can be seen in outline on the 
edge of the mass. Individual pyrite crystals infill pore spaces. 
Sample 80-9-5; Facies C; 442.5 m (1451.7 ft) 
Magnification: 156.25 X 


3. | Subangular metamorphic quartz grain with numerous sillimanite(?) needle-like 
inclusions. Oil stained clay lines the pore space and is found at grain-to-grain 
contacts. Small clay particles are seen within the pore space as well. Some may be in 
place, but most have been displaced during thin section preparation. 

Sample 80-9-1; Facies B; 439.5 m (1442 ft) 
Magnification: 390.63 X 
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Appendix |: Geophysical Well Logs for well 10-17-66-5W4 
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Appendix III: Calculation of Parameters Generated from Geophysical Well Log Data 
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The values obtained from the digitization of the geophysical 
well logs were used to generated parameters for shaly sand analysis. 


The equations used to obtain these values are as follows: 


(1) MestApresa 
st x 0.01 
b 
(2) N= Oy - Oye 
Fo GR 
(3) 6 A, Ana X Cp 
Ant ~ Ama 
(4) Gamma Ray Index = GR - GRoy 
GR», - GRoy 
(5) Bulk Volume Clay Fraction = 9. - 9, 
(6) Calculated Effective Porosity = (1 - q) 0 


C 


(7) q= 0. - 
05 


(8) Bulk Volume Oil = Winer Ol i " 
wt. %2 oil + wt. % water c 


(9) Bulk Volume Water = ( wt. % water )-o 
wt. 2 oil + wt. % water C 
where: A,¢ = sonic travel time for fluid (assumed 189 microsec/ft) 
A, = sonic travel time for point in question 


P, = bulk density for point in question 
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P. = bulk density of the fluid (assumed 1.0 gm/cc) 

by = neutron porosity of point in question 

= neutron porosity of fluid (assumed 100.0 percent) 
d, = corrected sonic porosity of point in question 


A. = sonic travel time for matrix: clean quartz sand 
(assumed 55 microsec/ft) 


Auf = sonic travel time for mud filtrate 
(assumed 189 microsec/ft) 
C = Compaction Correction: sonic travel time of formation shale _ 
100 
10as 
Too = Led 


GR = gamma ray value for point in question 
GRY = gamma ray value for clean quartz sand (assumed 20 API units) 


GR_, = gamma ray value for shale (132 API units) 


sh 
0. = core porosity 
0, = effective porosity 

q = volume of dispersed shale 


bo ="density porosity of point in question 
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